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19.  (cont.) 

The  NORSAR  Detection  Processing  system  has  been  operated  through  the 
period  1  April  -  30  September  1985  with  an  average  uptime  of  98.2 
per  cent.  Most  of  the  downtime  has  been  caused  by  disk  failure  at  the 
central  recording  computer.  A  total  of  2133  seismic  events  have  been 
reported  by  NORSAR  In  the  period.  The  performance  of  the  continuous 
alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC 
have  been  satisfactory.  After  a  redesign  of  the  system,  all  detection 
and  event  processing  Is  now  conducted  on  one  single  IBM  4341  computer, 
which  has  allowed  for  more  efficient  and  flexible  operation. 

Using  a  data  base  of  about  170  NORSAR-recorded  seismic  events  at 
regional  distances  during  the  past  15  years,  a  detailed  map  of  propa¬ 
gation  characteristics  of  short  period  regional  seismic  phases  has 
been  compiled.  Lg  shows  efficient  propagation  out  to  about  20  degrees 
at  1  Hz  and  10  degrees  at  4  Hz  over  continental  shield  paths  (e.g.,  to 
the  east  of  the  array),  whereas  Sn  Is  the  dominant  secondary  phase  for 
paths  that  cross  significant  tectonic  features  like  graben  structures 
and  continental  margins. (The  relative  amplitudes  of  P,  Sn  and  Lg  are 
strongly  frequency  dependent;  thus  Lg  Is  relatively  strong  at  low  fre¬ 
quencies  (around  1  Hz),  whereas  both  P  and  Sn  have  higher  dominant 
frequencies. 

-The  NORESS  on-line  detection  and  location  performance  have  been  eva¬ 
luated  both  for  regional  and  teleselsmlc  events.  A  case  study  of  132 
Western  Norway/North  Sea  events  (mainly  In  the  distance  range  300- 
600  km)  has  Indicated  a  detection  threshold  of  Ml  1. 5-2.0  for  NORESS. 
The  NORESS  eplcentral  solutions  are  In  most  cases  consistent  with 
those  provided  by  the  Western  Norway  seismic  network.  For  one  site 
where  a  large  number  of  known  explosions  (36)  were  processed,  the 
average  NORESS  location  error  was  10  km.  Occasional  larger  errors  are 
due  to  problems  in  the  automatic  phase  determination  at  NORESS  (e.g., 
Pg  Identified  as  Pn,  Sn  identified  as  Lg). 

A  statistical  study  of  NORESS  detections  during  April-September  1985 
(excluding  time  periods  with  transmission  line  problems)  has  shown  an 
average  of  about  100  detections  per  day,  of  which  about  20  per  cent 
are  teleselsmlc  P  phases,  about  25  per  cent  regional  P  phases  and  the 
remaining  55  per  cent  diave  slow  phase  velocities  (<  6  km/s)  corre¬ 
sponding  to  secondary  phases  or  possibly  noise  detections  .‘"'Most  of  the 
regional  phases  are  detected  during  local  day  time,  and  correspond  to 
quarry  blasts  and  other  local  explosions.  The  majority  of  P  phases  are 
detected  using  conventional  beams,  whereas  incoherent  (envelope)  beams 
are  very  efficient  in  detecting  secondary  phases. 

A  comparison  of  NORESS  and  NORSAR  has  shown  that,  for  the  six-month 
period  April-September  1985,  NORESS  detected  84.5  per  cent  of  the 
teleselsmlc  events  reported  in  the  NORSAR  monthly  bulletin  (only  days 
when  both  systems  were  operational  were  considered).  The  excellent 
teleselsmlc  capability  of  NORESS  is  further  confirmed  by  many  P-phases 
detected  at  NORESS  but  not  reported  by  NORSAR.  NORESS  estimates  of 
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19.  (cont.) 

phase  velocities  and  azimuths  are,  however,  not  as  accurate  as  those 
provided  by  NORSAR.  The  median  differences  In  NORSAR-NORESS  estimated 
azimuth  Is  12  degrees  (based  on  1376  common  events),  which,  based  on 
comparisons  with  PDE  solutions.  Is  mainly  attributed  to  uncertainty  in 
the  NORESS  estimates. 

A  study  of  P-wave  spectra  recorded  at  NORESS  has  been  initiated.  For 
Intermediate  and  low  magnitude  Eurasian  explosions  signal  energy  well 
above  the  noise  level  can  be  observed  for  frequencies  up  to  about 
10  Hz,  even  at  teleselsmlc  distances.  For  detection  purposes,  the  SNR 
(at  the  single-sensor  level)  Is  usually  best  in  the  2-4  Hz  band,  and 
the  beam  SNR  is  further  improved  by  selecting  appropriate  sub- 
geometries  of  NORESS  for  beamforming. 

An  efficient  method  for  the  calculation  of  ellipticity  corrections  has 
been  developed,  and  the  calculations  compare  well  to  results  published 
elsewhere.  Attenuation  of  seismic  energy  from  local  events  in  southern 
Norway  loas  been  studied  using  the  coda  decay  method  as  well  as 
spectral  ratios.  The  Q  values  are  found  to  be  comparable  to  values 
published  for  the  Canadian  Shield  at  high  frequencies  (12-15  Hz),  but 
are  more  strongly  frequency  dependent  than  Canadian  Shield  values. 

Various  signal  analysis  schemes  (three-component  analysis,  f-k  analy¬ 
sis  and  semblance)  have  been  tested  on  NORESS  data  from  a  presumed 
explosion  in  Western  Russia  and  have  been  found  to  produce  consistent 
and  reasonable  results  for  the  P  wavetraln.  The  f-k  method  has  the 
advantage  of  providing  efficient  scanning  of  the  whole  wavefleld 
space.  The  semblance  technique  is  very  robust,  but  requires  a  priori 
information  on  azimuth  for  easy  use  in  an  automated  mode.  3-component 
results  are  also  rather  stable  when  based  on  beam  traces  from  4 
stations. 

A  study  of  P-wave  coda  scattering  using  NORSAR  and  NORESS  data  has 
indicated  sources  of  Pg  scattering  mainly  to  the  north  and  east,  with 
a  clustering  that  appears  to  coincide  roughly  with  major  tectonic 
boundaries.  For  Pn  similar  results  are  found,  but  the  scattering 
sources  appear  to  be  more  dispersed  in  this  case.  _ 
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I 


SUMMARY 


This  Final  Technical  Report  describes  the  operation,  maintenance  and 
research  activities  at  the  Norwegian  Seismic  Array  (NORSAR)  for  the  period 
1  April  -  30  September  1985.  It  supplements  the  three  Semiannual  Technical 
Summaries  previously  submitted  during  the  contract  period,  which  has 
covered  1  October  1983  -  30  September  1985. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has  averaged 
98.2  per  cent,  during  April-September  1985,  which  is  the  same  as  for  the 
previous  half-year  period.  Disk  problems  in  August  contributed  strongly  to 
a  downtime  of  4.4%.  The  array  communications  system  for  02C  has  been  fre¬ 
quently  affected  during  the  period,  and  subarrays  02B,  03C  and  04C  also 
occasionally  •showed  high  error  figures.  Subarray  06C  resumed  operation  in 
July  and  now  has  an  acceptable  performance.  A  total  of  2133  events  were 
reported  in  this  period,  giving  a  daily  average  of  11.7  events. 

During  this  period  the  system  design  has  been  changed,  so  that  the  entire 
system  including  data  acquisition,  detection  processing,  event  processing 
and  plotting  is  now  performed  by  only  one  IBM  4341  computer.  This  change 
will  facilitate  both  operation  and  maintenance  in  the  future. 

Much  of  the  field  staff's  work  during  April  and  May  concerned  preparations 
for  the  official  NORESS  opening,  which  took  place  in  June.  During  the 
months  of  July-September  the  field  maintenance  center  was  moved  from  Stange 
to  Hamar.  In  the  reporting  period,  field  maintenance  work  was  divided 
fairly  evenly  between  NORSAR  and  NORESS. 

The  research  activity  is  summarized  in  section  VII.  Section  Vll.l 
discusses  propagation  characteristics  of  regional  phases  recorded  at 
NORSAR.  Section  VII. 2  contains  an  evaluation  of  NORESS  real  time  processing 
performance  with  a  discussion  of  a  case  study  for  132  Western  Norway/ 

North  Sea  events.  Section  VII. 3  presents  the  initial  results  from  NORESS 
detection  processing,  while  Section  VII. 4  contains  a  comparison  of  the 
NORESS  and  NORSAR  processing  systems.  Section  VII. 5  presents  a  study  of  P- 
wave  spectra  from  NORESS  recordings.  Section  VII. 6  presents  a  computational 


method  for  the  calculation  of  ellipticity  corrections.  Section  VII. 7  pre¬ 
sents  a  study  of  the  attenuation  of  seismic  energy  from  local  events  in 
southern  Norway.  In  section  VII. 8  a  new  technique  for  3-component 
seismogram  analysis  is  presented,  while  Section  VII. 9  discusses  wavefield 
decomposition  schemes  based  on  NORESS  regional  event  records.  Section 
VII. 10  contains  a  discussion  of  the  on-going  efforts  to  determine  the 
nature  of  the  teleseismic  P  coda  seen  on  short  period  seismograms. 


There  have  been  92  breaks  in  the  otherwise  continuous  operation 
of  the  NORSAK  online  system  within  the  current  b-month  reporting 
interval.  The  uptime  percentage  for  the  period  is  98.2,  which  is 
the  same  as  for  the  previous  half-year  period. 

Fig.  II. 1.1  and  the  accompanying  Table  11.1.1  both  show  the  daily 
DP  downtime  for  the  days  between  1  April  and  30  September  1985. 
The  monthly  recording  times  and  percentages  are  given  in  Table 
II. 1.2. 


The  breaks  are  grouped  as  follows; 

a)  Computer  hardware  failure  11 

b)  Stops  related  to  programming  work 

or  error  6 

c)  Hardware  maintenance  stops  27 

d)  Power  merges  and  power  breaks  6 

e)  TOD  correction  29 

f)  Communication  line  failures  13 


The  total  downtime  for  the  period  was  77  hours  and  44  minutes. 

The  mean-time-between-f ailure  (MTBF)  was  2.0  days,  as  compared  to 
1.4  for  the  previous  period. 


J.  Torstveit 


Detection  Processor  downtime  in  the  period  1  April  -  30  September  1985. 
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Daily  DP  downtime  in  the  period  1  April  -  30  September  1985. 
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Month 


DP  uptime  DP  uptime 
hours  % 


No.  of 
DP  breaks 


APR 

700.40 

97.3 

24 

MAY 

743.15 

99.9 

10 

JUN 

712.67 

99.0 

14 

JUL 

735.05 

98.8 

14 

AUG 

711.42 

95.6 

14 

SEP 

711.58 

98.8 

92 

No.  of  days  DP  MTBF* 
with  breaks _ (days) 
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Table  II. 2.1  reflects  the  performance  of  the  system  throughout 
the  reporting  period.  With  reference  to  the  table  we  can  see  that 
02C  has  been  frequently  affected.  02B,  03C  and  04C  also  occa¬ 
sionally  showed  high  error  figures. 


06C  resumed  operation  12  July  after  having  been  connected  to  the 
PCM  system  between  Rokosj0en  telephone  central  and  Hamar,  and  now 
has  an  acceptable  performance. 


Most  the  irregularities  have  been  caused  by  faulty  and/or  damaged 
cables  C02B,  02C,  03C,  04C),  changes  and/or  expansion/reconfigura¬ 
tion  at  telestations  (Lillehammer/02C) ,  lightning  (especially 
exposed  04C,  06C),  Modcomp  problems  during  resynchronizing 
sequences  (especially  when  in  connection  with  deteriorated 
line(s)),  and  power  outages  (04C). 


Summary 

Apr  85  02C  affected  in  connection  with  reconfiguration/ 

expansion  at  NTA/Lillehammer  (between  1  and  12  Apr) 

04c  was  down  due  to  a  bad  communication  line  (9-12  Apr). 

06C  comm,  system  still  experiencing  problems. 

In  a  letter  to  NTA,  NORSAR  recommended  reconnection 
to  original  line  between  Hamar  and  06C. 


May  85  This  was  a  relatively  good  period  in  spite  of  02C  outa¬ 
ges  of  25%. 

An  attempt  to  strap  the  06C  NOKIA  modems  to  external 
synchronization  did  not  improve  the  performance  of  that 


subarray.  In  internal  synchronization,  spikes  created 
unacceptable  conditions. 

Jun  85  Subarrays  OlA,  OIB,  02B  and  03C  nearly  error  free  as 
they  were  in  the  previous  two  months. 

02C  was  down  3, 4-6, 7  and  8-10  June,  last  outage  12  June. 

04c  was  affected  between  3  and  4  June  as  we  were  not 
able  to  resynchronize  the  system. 

06C  was  still  out  of  operation,  awaiting  reconnection 
to  the  original  line. 

Jul  85  Only  three  subarrays  were  nearly  free  from  errors: 
01A-02B. 

02C-04C  had  high  error  figures. 

06C  was  reconnected  to  the  original  line  12  July,  after 
which  remarkable  improvement  was  observed. 

The  error  source  which  caused  the  varying  02C  perfor¬ 
mance  was  not  stated  by  NTA  this  period;  a  bad  connec¬ 
tion  or  an  intermittently  operating  item  may  have  caused 
the  frequent  outages.  Week  29  NTA  technicians  visited 
the  CTV  and  were  able  to  reinstate  normal  operation. 

A  faulty  cable  between  Osmoen  and  the  subarray  03C 
caused  reduced  performance  (week  31). 

Lightning  caused  power  problem  (week  28)  at  04C. 

Aug  85  In  August  NTA/Hamar  started  to  repair  a  cable  near 

Mesnalien,  which  also  contained  the  02C  comm,  cable. 
Almost  simultaneously  the  same  cable  was  damaged  two 
other  places,  one  was  located  near  the  CTV.  This 
situation  created,  of  course,  special  problems  and 
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delayed  the  termination  of  the  splicing  job. 

03C  was  down  two-thirds  of  week  32,  due  to  a  cable  of 
reduced  quality;  part  of  it  was  replaced. 

02B  was  affected  week  35,  probably  caused  by  a  mistake 
by  NT A/ Hama r. 

Lightning  caused  OAC  isolation  transformer  damage  13 
August,  and  consequently  loss  of  power.  Tranformer  was 
replaced  21  August  and  the  subarray  resumed  operation. 

Besides,  synchronization  problems  affected  the  comm, 
system  performance  for  02C-06C. 

In  order  to  establish  normal  conditions  it  was  necessary 
to  stop  and  restart  the  communication  machine  (Modcomp ) 
and/or  manipulate  isolation  switches  between  the  modems 
and  the  comm,  machine,  specially  week  35.  Such  problems 
may  arise  in  cases  where  the  comm,  machine  is  heavily 
loaded  and  tries  to  synchronize  a  comm,  line  with 
reduced  quality. 

Sep  85  Communication  systems  specially  affected  in  September 
were  02B,  03C  and  U4C. 

Subarray  02B  lost  power  29  September  and  remained 
inoperatlonal  the  rest  of  the  month. 

03C  was  down  between  2  and  6  September.  In  the  same 
period  04C  caused  heavy  Modcomp  load  (repeated 
resynchronization  attempts)  due  to  a  cable  fault  near 
the  CTV  causing  gradually  reduced  performance. 

06C  was  also  affected  by  synchronization  problems  (2 
September) . 


O.A.  Hansen 
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IV.  IMPROVEMENTS  AND  MODIFICATIONS 

NORSAR  on-line  system  using  4341  and  MODCOMP  Classic 


We  refer  to  the  detection  processor  operation  statistics  for 
detailed  information  about  uptime,  communication  and  processing 
performance. 

Processing  system: 

The  NORSAR  online  system  has  previously  been  using  two  computers  to  do 
detection  processing,  generate  event  bulletins  and  plot  events.  The 
system  design  has  now  been  changed  so  that  everything  from  data 
acquisition  to  event  plotting  is  achieved  using  only  one  IBM  4341  com¬ 
puter.  It  is  easier  for  us  to  both  maintain  and  operate  the  system 
when  only  one  computer  is  involved. 

The  following  functions  are  implemented  on  the  NORSAR  on-line 
computer: 


Functions  Program 

Communication  with  MODCOMP  (DPCOM) 

Detection  processing  (DPDAQ) 

and  storing  of  data 

Event  processing,  plotting  (DPEP) 

Automatic  preparation  of  bulletins  (AFTAC) 


The  computer  used  for  these  functions  is  an  IBM  Model  4341  LOl. 

The  average  system  load  is  about  28%.  Thus,  there  is  still  computing 
power  left  even  with  the  current  full  operation  for  implementing  new 
processing  algorithms. 


Data  fron  06C 


The  06C  subarray,  which  for  several  years  was  used  in  experimental 
configurations  related  to  NORESS  developments,  has  now  been  reestab¬ 
lished  in  its  original  configuration,  as  part  of  the  standard  NORSAK 
array. 

Data  quality 

From  time  to  time  we  have  communication  line  problems  which  give  us 
a  significant  number  of  false  detections.  In  the  automatic  bulletin 
transferred  to  AFTAC  by  telex,  the  number  of  DPX  entries  during  such 
periods  is  very  large.  Also  the  automatic  event  lists  (EPX)  are 
affected,  since  the  current  processing  system  is  not  able  to  mask  all 
such  occurrences.  We  are  working  on  improvements  in  regard  to  this 
problem. 

ARPANET 

The  upgrade  of  our  ARPANET  connection  is  planned  for  completion 
during  November  85.  Our  experience  using  ARPANET  for  NORESS  bulletin 
transfer  to  the  CSS  is  very  good  and  will  be  even  more  flexible  when 
the  planned  ARPANET  host  connection  is  established. 

NORESS 

A  significant  part  of  our  research  under  the  AFTAC  contract  is  con¬ 
nected  to  studies  of  regional  seismic  phases,  using  data  from  the  new 
NORESS  array. 


V. 


MAINTENANCE  ACTIVITIES 


V « 1 _ Activities  In  the  field  and  at  the  Maintenance  Center 

Visits  to  the  NORSAR  array  were  in  May  and  June  limited  to  actions 
with  high  priority,  as  the  preparations  for  the  official  opening  of 
the  NORESS  array  3  June  needed  all  resources  available  at  the 
NMC. 

In  July,  August  and  September  normal  maintenance  work,  was  resumed. 

A  great  deal  of  time  was  also  spent  on  planning/preparations 
in  connection  with  moving  to  the  new  maintenance  center  in  Hamar. 

SP  and  LP  instruments  have  routinely  been  remotely  monitored  from 
NDPC  and  LP  seismometers  adjusted  whenever  needed. 

A  summary  of  the  activities  performed  during  the  period  is  given 
in  Table  V.l.  For  completeness,  we  list  activities  carried  out  at 

both  the  NORSAR  and  NORESS  sites. 

Subarray/ 

area _ Task _ Date 

NORESS  A  satellite  transmitter  amplifier  (TWT)  24  April 

failed  and  was  replaced  immediately. 

02C  An  attempt  to  replace  an  RA-5  SPOl  failed  25  April 

due  to  bad  weather  conditions. 


Task 


Date 


Subarray/ 

area 


NORESS/ 

NMC 

In  May  and  June  most  time  was  spent  on  pre¬ 
parations  and  planning  in  connection  with  the 
NORESS  official  opening  3  June  and  the  move 
to  the  new  NORSAR  maintenance  center  in 

Hamar.  Also  in  June  the  NMC  staff  assisted 
Sandia  representatives  in  modifying  equipment 
at  the  earth  station  (4-27). 

May/ June 

NORESS 

10  channels  in  the  NORESS  array  were  out  of 
operation  a  few  hours  due  to  a  short  circuit 
in  a  junction  box  which  was  repaired  by  the 

NMC  staff. 

8  July 

NORESS 

NORESS  earth  station  was  down  approx.  6  hours  15  July 
due  to  a  high  voltage  transformer  damaged  by 
lightning.  The  NMC  staff  was  engaged  in  getting 
the  station  operational  again. 

NORESS 

Lightning  caused  power  trouble  and  NMC 
engagement.  The  station  was  down  for  about 

2  hours  and  15  minutes. 

23  July 

NORESS 

In  connection  with  installation  of  new 
switches  on  High  Voltage  Transformer  (31  July) 
the  station  was  down  for  about  2  hours.  NMC 
took  part  in  the  work. 

31  July 

06C 

NMC  staff  collected  former  06C  experimental 
vaults  and  brought  them  to  the  center  station. 

July 

OIB 

Visited  and  replaced  RA-5  amplifier  SP02,05. 
Alto  localized  faulty  cable  (common  for  SPOl, 
05),  and  in  addition  a  broken  cable  to  SP02. 
Both  cable  faults  caused  by  lightning. 

5  August 

OIB 

Repairs  and  cable  splicing. 

19,20,21 
&  26  August 

02C 

Replaced  both  Remote  Centering  Devices  on 

EW  seismometer  (LP). 

Bad  comm,  line  localized. 

13  August 

14  August 

0  3C 

Replaced  FP  Remote  Centering  Device  on 

Vertical  and  NS  seismometers  (LP). 

14  August 
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Subarray/ 

area 

Task 

Date 

04C 

Power  isolation  transformer  (NORSAR  property) 
was  damaged  by  lightning.  NMC  spent  three 
days  on  this  job  (incl.  demounting  ex.  05C) 
and  transport. 

August 

NORESS 

NORESS  central  station  was  visited  in 
connection  with  analyses  of  the  Ay  HF-channel 
spectrum. 

August 

02B 

(Telem. ) 

All  02B  (Telemetry)  stations  have  been 
visited,  batteries  and  seismometer  cables 
have  been  replaced. 

September 

NMC 

Engagements  in  September  have  been  concen¬ 
trated  on  preparing  equipment  for  a  new 
experimental  small  array. 

September 

Table  V.i  Activities  in  the  field  and  at  the  NORSAK  Maintenance 
Center,  April  -  September  1985. 


V.2 _ Improvements  and  modifications 


i 


No  NORSAR  equipment  modification  was  carried  out  during  the  reporting 
period. 


As  of  30  September  1985  the  following  channels  deviated  from 
tolerances : 


04C 

01, 

02,  03,  06  and  07 

06C 

03, 

08 

OlA 

01 

8  Hz  filter 

02 

8  Hz  filter 

04 

Attenuated 

Note  that  the  deviations  for  OlA  do  not  reflect  errors,  as  these 
channels  are  used  experimentally.  The  deviations  at  OAC  and  06C 
will  be  corrected. 


U.A.  Hansen 
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Wessel,  P.  &  E.S.  Husebye  (1985):  The  Oslo  graben  gravity  high 
and  taphrogeneals,  submitted  for  publication 


L.B.  Loughran 


at  NORSAR 


About  170  regional  events  recorded  at  NORSAR  during  the  period 
1971-1985  have  been  analyzed  with  respect  to  propagation  charac¬ 
teristics  for  the  phases  Sn  and  Lg  in  particular.  The  motivation 
for  our  study  was  to  obtain  calibration  data  for  the  real  time 
processing  of  data  recorded  on  the  new  NORESS  array.  In  addition, 
it  is  expected  that  such  a  study  will  shed  some  light  on  the 
general  relationship  between  characteristics  of  regional  seismic 
phases  and  large-scale  geological  features. 

The  event  location  algorithm  of  the  current  version  of  the  NORESS 
online  processing  package  (RONAPP)  is  based  on  the  assumption 
that  the  largest  secondary  phase  detected  is  Lg,  propagating  at  a 
group  velocity  of  3.5  km/s.  This  assumption  is  valid  for  many 
propagation  paths  to  regional  distance  from  the  NORSAR  and  NORESS 
arrays.  There  are,  however,  notable  exceptions  to  this  simplified 
picture,  as  illustrated  in  Fig.  VII. 1.1. 

Data  analysis  procedure 

For  the  170  events  in  the  data  base,  event  files  of  appropriate 
lengths  to  include  all  phases  of  interest  were  copied  from  the 
tape  archive.  In  order  to  minimize  possible  effects  from  the 
geology  near  the  receiver  site,  the  analysis  was  carried  out  on 
data  from  the  same  seismic  channel  (02B01)  for  nearly  all  events. 
Fig.  VII. 1.2  shows  an  example  of  the  output  resulting  from  the 
processing  of  events  in  our  data  base.  The  seismograms  in  this 
figure  are  shown  for  unflltered  data  (top)  and  three  band-pass 
filters.  We  needed  criteria  for  picking  phase  onsets,  which  are 
Independent  of  the  human  analyst's  personal  judgement.  To  this 
end,  we  computed  and  plotted  the  logarithm  of  the  STA/LTA-ratio 


for  each  trace.  The  STA/LTA-algorithm  used  here  for  a  single  sen¬ 
sor  Is  equivalent  to  the  one  used  In  the  detection  procedure  In 
the  RONAPP  processing  package,  where  STA/LTA  Is  computed  for  a 
number  of  beams.  In  this  manner,  our  analysis  ties  In  with  the 
current  automatic  processing  of  NORESS  data.  For  this  study,  we 
specified  an  STA/LTA  threshold  of  3.0  for  acceptance  of  a  signal 
as  a  true  phase  arrival.  This  level  would  correspond  to  a 
somewhat  higher  (and  operationally  acceptable)  detection 
threshold  for  the  real  time  array  processing,  due  to  the  SNR  gain 
offered  by  the  array  beamforming. 

Now,  our  analysis  procedure  has  been  to  go  through  170  plots  of 
the  kind  shown  in  Fig.  VII. 1.2  and  characterize  the  arrivals  in 
the  S  wavetrain  for  each  of  the  three  filter  bands.  The  case  of 
dominant  Lg  (over  Sn)  was  subdivided  into  three  different  codes 
as  follows: 

(1)  Only  Lg  can  be  seen  in  the  S  wavetrain. 

(2)  Lg  Is  the  dominating  phase,  but  Sn  can  be  discerned. 

(3)  Lg  is  still  the  dominating  phase,  but  Sn  is  not  very 
much  smaller. 

An  equivalent  subdivision,  also  on  a  three-point  scale,  was 
adopted  for  the  case  of  dominant  Sn.  Finally,  a  code  was  assigned 
to  the  case  of  comparable  Lg  and  Sn  amplitudes.  Phases  with  cor¬ 
responding  group  velocities  in  the  interval  3.80-4.70  km/s  were 
classified  as  Sn,  with  Lg  corresponding  to  group  velocities  In 
the  range  3.30-3.80  km/s. 

Analysis  results  —  source  area  aapping 

Fig.  VII. 1.3  shows  the  results  derived  from  analysis  of  the  data 
In  our  data  base.  In  this  figure,  we  assign  the  character  of  the 
S  wavetrain  as  observed  In  the  NORSAR  seismometer  record,  to  the 


source  area.  We  have  confined  ourselves  to  the  use  of  three  codes 
only  in  this  figure  —  corresponding  to  (1)  dominant  Lg,  (2)  com¬ 
parable  Lg  and  Sn  and  (3)  dominant  Sn.  Areas  with  no  events  or 
for  which  the  available  data  were  inconclusive  with  respect  to 
appropriate  characterization  of  the  S  wavetrain  appear  unshaded 
in  the  figure.  Color-coded  maps  including  all  7  codes  defined 
above  have  also  been  prepared.  These  maps  offer  more  details  and 
often  show  gradual  transitions  from  dominating  Lg  to  dominating 
Sn  in  certain  directions.  Several  conclusions  can  be  drawn  from 
Fig.  VII. 1.3: 

For  an  area  covering  almost  all  of  Norway,  Sweden, 
Finland  and  the  Baltic  Sea,  Lg  is  consistently  the 
largest  secondary  phase  for  all  three  frequency  bands. 

1  Hz  Lg-waves  are  clearly  observed  at  NORSAR  for  sources 
in  all  of  Poland  and  major  parts  of  the  western  U.S.S.R. 
up  to  distances  of  about  20°. 

Sn  is  the  dominant  phase  for  events  from  the  U.K., 
Belgium,  Iceland  and  the  North  Atlantic  mid-ocean  ridge, 
and  for  events  located  to  the  west  of  the  escarpment  of 
the  Norwegian  continental  shelf. 

The  maps  tend  to  be  richer  in  details  for  higher  fre¬ 
quencies,  with  more  pronounced  boundaries  between  source 
regions . 

Wave  propagation  characteristics  related  to  tectonic  features 

The  map  of  Fig.  VII. 1.4  shows  the  tectonic  features  that  ten¬ 
tatively  could  be  related  to  the  characteristics  of  propagation 
of  regional  phases  in  the  area  surrounding  the  NORSAR  array. 
Indeed,  we  find  a  very  close  relationship  between  the  tectonic 


features  and  the  character  of  the  S  wavetrain.  By  comparing  the 
maps  in  Figs.  VII. 1.3  and  VII. 1.4  and  bearing  In  mind  also  the 
significance  of  unshaded  areas  in  Fig.  VII. 1.3,  the  following 
points  can  be  made: 

-  Lg  wave  propagation  is  very  efficient  for  propagation 

paths  that  are  confined  to  continental  shield  type 
structures  and  not  intersected  by  major  tectonic  units. 
In  these  directions,  1  Hz  Lg  waves  propagate  to  20°  from 
the  source,  while  at  4  Hz  Lg  is  seen  out  to  about  10° 
from  NORSAR. 

No  Lg  is  observed,  or  Sn  dominates  over  Lg  for  propaga¬ 
tion  paths  that  cross  significant  tectonic  features  like 
graben  structures  and  continental  margins. 

Particularly  at  4  Hz  there  is  a  close  connection  between 
some  of  the  finer  tectonic  details  and  the  propagation 
characteristics.  Since  Lg  is  a  phase  mostly  composed  of 
waves  trapped  in  the  crust,  this  phase  is  likely  to  be 
highly  susceptible  to  changes  in  the  crustal  structure. 
This  effect  should  be  more  pronounced  at  higher  frequen¬ 
cies,  i.e.,  for  shorter  wavelengths.  At  4  Hz,  even  a 
feature  like  the  Tornquist  line  seems  to  be  an  effective 
barrier  to  the  propagation  of  Lg. 

Implications  for  real  time  data  processing 

The  basic  assumption  in  RONAPP  that  for  regional  events,  Lg  is 
the  strongest  phase  in  the  S  wavetrain,  is  not  valid  throughout 
the  area  of  Interest.  However,  the  maps  in  Fig.  VII. 1.3  will 
enable  the  analyst  to  assess  the  likelihood  of  a  detected  secon¬ 
dary  phase  being  either  Sn  og  Lg,  when  reviewing  the  results  from 
the  automatic  event  processing.  On  this  basis,  the  eplcentral 
distance  estimate  could  be  corrected,  if  necessary.  In  some 


cases,  the  phase  identification  will  still  be  ambiguous,  and 
research  is  continuing  towards  developing  further  criteria. 

A  next  step  would  be  to  utilize  the  regionalization  contained  in 
this  contribution  for  an  improvement  of  the  real  time  processing. 
We  think  that  the  results  from  this  study  and  knowlegde  about  the 
general  seismicity  pattern  of  the  area  together  with  the  infor¬ 
mation  already  available  from  the  detection  processing  (direction 
of  arrival,  signal  frequency  and  possibly  travel  time  differences 
between  several  arrivals  in  the  S  wavetraln),  constitute  a  basis 
for  such  improvements.  Eventually,  this  information  may  assist  in 
an  expert  system's  approach  to  regional  array  processing. 

T.  Kvaerna 

S.  Mykkeltveit 
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Illustration  of  variation  of  relative  importance  of 
the  phases  Sn  and  Lg.  The  standard  group  velocities 
of  4.5  and  3.5  km/s,  commonly  assigned  to  Sn  and 
Lg,  respectively,  are  marked  by  dashed  lines.  The 
upper  three  traces  cover  the  distance  interval 
480-550  km,  while  the  lower  three  traces  correspond 
to  epicentral  distances  in  the  range  1225-1320  km. 
The  location  of  the  NORSAR  array  is  denoted  by  a 
ring  on  the  map,  and  the  traces  are  from  NORSAR 
seismometer  02B01.  The  data  are  bandpass  filtered 
1-5  Hz. 


Fig.  VII. 1.3  Characterization  of  Lg  versus  Sn  propagation  effi¬ 
ciency,  for  the  three  different  filter  bands.  The 
shading  codes  are  as  follows:  (1)  Lg  dominant  phas 
in  the  S  wavetrain;  (2)  Lg  and  Sn  comparable  and 
(3)  Sn  dominant.  Note  that  the  codes  are  assigned 
to  the  source  area,  and  all  characterizations  are 
related  to  seismograms  as  they  appear  on  NORSAR 
(filled  circle)  seismometers. 


Fig.  VII. 1.4  Map  showing  geological  features  that  are  related  to 
characteristics  of  propagation  of  regional  phases 
to  the  NORSAR  array.  Shading  codes  are  (1)  Oceanic 
or  oceanic  type  crust  in  the  North  Atlantic;  (2) 

The  Iceland-Faeroe  Ridge;  (3)  Graben  structures; 

(4)  Major  mountain  belts;  and  (5)  The  Tornquist 
line.  Unshaded  areas  (except  in  the  Mediterranean 
region)  are  of  continental  type  without  significant 
orogenic  features. 
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VII. 2 _ Evaluation  of  NQRESS  real  time  processing  performance: 

Case  study  for  132  Western  Norway/North  Sea  events 

The  performance  of  the  real  time  processing  of  regional  events 
detected  and  located  by  the  NORESS  array  is  currently  being 
evaluated.  For  this  evaluation  it  is  crucial  to  be  able  to  com¬ 
pile  a  data  base  of  events  that  have  been  reliably  located  by 
another  agency.  Specifically,  one  would  look,  for  a  region  with  a 
fair  number  of  seismic  events  and  also  well  covered  by  seismic 
stations.  One  such  region  within  regional  range  from  NORESS  is 
Western  Norway. 

During  the  period  April-September  1985  a  total  of  132  seismic 
events  were  detected  and  located  by  the  RONAEP  real  time  pro¬ 
cessing  package  for  NORESS  data,  and  also  reported  by  the  Western 
Norway  Seismic  Network.  Fig.  VII. 2.1  shows  the  location  of  the 
NORESS  array,  the  location  of  the  six  stations  constituting  the 
network,  and  the  area  (hatched)  within  which  the  seismic  events 
occurred.  These  events  are  typically  presumed  earthquakes 
offshore  and  presumed  explosions  onshore.  The  magnitudes  are  in 
the  range  Ml  1.5-3. 5,  and  the  epicentral  distances  relative  to 
NORESS  range  from  200  to  nearly  700  km.  41  of  the  events  were 
verified  through  reports  as  explosions  and  their  locations  are 
known.  These  explosions  occurred  at  site  1  (3  events),  site  2 
(36  events)  and  site  3  (2  events).  There  is  a  concentration  of 
events  in  the  coastal  area  at  around  60°N,  in  addition  to  the 
explosion  site  2. 

Data  analysis 

The  real  time  processing  of  NORESS  data  produces  standardized 
plots  which  provide  the  analyst  with  all  essential  information  on 
the  performance  of  the  automatic  event  processor.  Fig.  Vll.2.2 
gives  an  example  of  such  plots  tor  one  regional  event.  The  loca¬ 
tion  of  a  regional  event  results  from  consecutive  detections  of 
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P-  and  S-type  arrivals  (identified  as  P  or  S  according  to  phase 
velocity)  with  a  common  azimuth.  The  range  is  estimated  from  the 
travel-time  difference  between  the  first  P  phase  (usually  Pn)  and 
the  strongest  phase  in  the  S  wavetrain,  which  is  assumed  to  be 
Lg,  propagating  at  a  constant  group  velocity  of  3.5  km/s.  RONAPP 
utilizes  the  bearing  estimated  for  the  Lg  phase  in  the  event 
location  procedure. 


Fig.  VII. 2. 3  shows  differences  bdtween  P  and  Lg  azimuths  for  all 
events  in  our  data  base.  The  median  of  the  absolute  values  of 
these  differences  is  5  degrees.  This  is  an  estimate  then  of  the 
consistency  that  should  be  expected  among  the  azimuths  automati¬ 
cally  determined  for  the  various  phases,  for  the  current  choice 
of  processing  parameters.  Our  P  and  Lg  azimuths  were  compared 
with  those  corresponding  to  the  locations  by  the  network.  We 
found  that  for  47  events  the  Lg  azimuth  came  closest  to  the 
"true"  value.  For  30  events  the  P  azimuth  was  the  better,  and  for 
the  remaining  55  events,  the  P  and  Lg  azimuth  deviations  from  the 
true  values  were  about  the  same. 


Fig.  VII. 2. 4  shows  the  differences  between  the  ranges  estimated 
by  the  NORESS  online  processing  and  the  ones  corresponding  to  the 
network  solutions.  As  can  be  seen  from  this  figure,  there  is  a 
tendency  by  the  NORESS  event  location  procedure  to  underestimate 
the  range  for  events  in  this  region.  In  addition  to  the  events 
included  in  Fig.  VII. 2. 4  there  are  35  events  in  the  data  base 
with  known  location  at  site  2  (Fig.  VII.2.1),  for  which  this 
range  difference  is  between  -10  and  -20  km. 

Bias  sources  and  location  improvements 

Both  systematic  and  random  errors  contribute  to  the  deviations  in 
event  locations  reported  above.  The  location  uncertainty  of  the 
network  must  be  considered  a  random  error  in  this  context.  The 


network's  average  location  errors  for  the  explosions  with  known 
epicenters  are  as  follows:  30  km  for  site  1,  10  km  for  site  2  and 
7  km  for  site  3.  For  a  number  of  the  events  in  our  data  base, 
rather  high  RMS  residuals  associated  with  the  network  solutions 
indicate  that  Lg  arrivals  may  have  been  mistaken  for  Sn  in 
reading  phase  arrivals.  This  warrants  a  network  relocation 
experiment  for  all  events  in  our  data  base,  taking  advantage  of 
the  regionalization  described  in  section  VII. 1  of  this  volume.  We 
think  that  such  a  relocation  experiment  could  result  in  improved 
network  locations. 

The  frequency-wavenumber  analysis  leading  to  an  estimate  of  the 
arrival  azimuth  is  performed  on  a  finite  grid  of  41  x  41  points 
in  wavenumber  space.  This  grid  size  limits  the  resolution  to 
approximately  3°  for  the  slower  phases  like  Lg  and  typically  5° 
for  regional  P  phases.  An  increase  in  the  number  of  grid  points 
or  a  restriction  of  the  frequency-wavenumber  analysis  to  a 
smaller  part  of  the  wavenumber  space  would  increase  the  resolu¬ 
tion  power  of  the  azimuth  estimator. 

The  travel  time  model  which  is  part  of  the  event  location 
algorithm  in  the  RONAPP  processing  package  is  based  on  a  uniform 
Moho  depth  of  40  km.  This  is  appropriate  for  a  majority  of 
regional  propagation  paths  to  NORESS.  For  the  events  being 
studied  here,  however,  a  crustal  thickness  of  30  km  in  the  source 
region  and  35  km  at  the  NORESS  array  would  be  adequate.  Our  loca¬ 
tion  procedure  assumes  surface  sources,  while  a  number  of  the 
events  in  our  data  base  are  earthquakes,  of  unknown  but  probably 
intermediate  crustal  depth.  With  these  changes  in  our  travel  time 
model,  the  automatically  determined  ranges  would  Increase  by 
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10  km  for  surface  sources  and  by  15-20  km  for  10-15  km  deep 
earthquakes . 

Careful  inspection  of  the  NORESS  seismograms  of  our  data  base 
resulted  in  the  following  observations: 

For  31  events  in  the  distance  range  350-500  km  we  found 
reason  to  believe  that  the  Pn  phase  had  gone  undetected. 
For  these  cases  the  P  phase  detected  is  mostly  likely 
Pg.  For  most  of  these  events,  Pn  can  be  discerned  as  a 
weak  arrival  preceding  the  much  stronger  Pg.  The  time 
difference  is  of  the  order  of  3.5-5  sec,  corresponding 
to  an  underestimation  of  epicentral  distance  of  the 
order  of  20-30  km.  For  28  out  of  these  31  Pg  arrivals, 
the  phase  velocity  is  in  the  range  6. 5-7. 5  km/s,  which 
is  considered  very  adequate  for  this  phase.  For  the 
three  remaining  signals,  the  phase  velocity  exceed 
8  km/s,  which  would  normally  be  indicative  of  a  Pn 
arrival.  This  result  indicates  a  potential  of  automati¬ 
cally  discriminating  between  Pn  and  Pg. 

For  9  events,  the  epicentral  range  was  underestimated  by 
more  than  70  km.  All  of  these  cases  are  related  to 
mistaking  Sn  for  Lg,  with  Sn  being  the  dominant  secon¬ 
dary  phase  or  the  only  secondary  phase  detected.  With 
proper  phase  assignments,  the  ranges  for  these  9  events 
would  have  been  correctly  determined. 

The  nvin  reason  tor  overestimating  the  range  is  late 
detection  ot  tlie  Lg  phase. 


Corrections  along  the  lines  indicated  would  lead  to  significant 
improvements  in  the  range  estimation  for  the  Western  Norway/North 
Sea  events. 

More  generally,  we  think  that  the  experience  gained  through  this 
case  study  will  result  in  substantial  improvements  in  the  perfor¬ 
mance  of  the  online  processing  of  regional  events  recorded  on 


NORESS. 


« 


Fig.  VII. 2.1 


Map  showing  the  location  of  the  NORESS  array  (encircled 
station  symbol),  the  sixt  stations  of  the  Western  Nor¬ 
way  Seismic  Network,  the  three  locations  ("1",  "2"  and 
"3")  for  reported  explosions,  and  the  source  region 
(hatched)  with  the  132  seismic  events  in  our  data  base. 
The  location  of  some  of  the  offshore  earthquakes  was 
constrained  also  by  signals  recorded  on  seismic  stations 
to  the  west  of  the  source  area,  i.e.,  on  Shetland  and 
in  Scotland. 
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Fig.  VII. 2. 2a  Individual  NORESS  traces  for  one  of  the  events  in  our 

data  base.  The  panel  covers  min  of  bandpass  filtered 
records  (3-5  Hz).  The  four  arrows  on  top  of  the  records 
indicate  detections  by  the  RONAPP  real  time  processing 
package.  The  broad  arrows  correspond  to  Pn  and  Lg  arri¬ 
vals,  which  have  been  combined  to  locate  the  event  at 
the  geographical  position  given  on  top  of  the  figure. 
The  two  additional  detections  correspond  to  Pg  and  Sn 
arrivals.  P  and  Lg  beams  are  shown  in  the  bottom  part 
of  the  figure. 
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Fig.  VII. 2. 2b  NORESS  automatic  analysis  results  for  the  event  of 

Fig.  VII. 2. 2a.  The  Pn  phase  (left)  and  Lg  phase  (right) 
are  shown  for  all  instruments.  Frequency-wavenumber 
spectra  are  shown  for  Pn  (top)  and  Lg  (bottom).  These 
spectra  are  computed  for  phase  identification  (P-  or  S- 
type  arrival  according  to  phase  velocity)  and  estima¬ 
tion  of  arrival  azimuth. 
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Fig.  VII. 2. 4  The  histogram  shows  the  differences  between  the  ranges 
estimated  by  the  NORESS  online  processing  and  the  ones 
corresponding  to  the  network  solutions,  for  the  events 
in  our  data  base.  35  events  with  known  location  at 

2  (Fig.  VII. 2.1),  for  which  this  difference  was 
between  -10  and  -20  km,  are  not  included  in  the  figure 
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VII «3 _ Initial  results  from  NORESS  detection  processing 

Since  January  1985  data  from  the  new  small-aperture  array  NORESS 
in  Norway  have  been  processed  in  real  time  at  the  NORSAR  data 
center  at  Kjeller.  The  data  used  in  the  detection  processing 
comprise  25  SPZ  channels,  deployed  over  an  area  3  km  in  aperture 
and  sampled  at  40  Hz  rate.  The  detection  algorithm  (RONAPP)  has 
been  described  by  Mykkeltveit  and  Bungum  (1984),  and  the  current 
implementation  consists  of 

Digital  narrow-band  filtering  using  a  set  of  Butterworth 
recursive  bandpass  filters 

Beamforming,  using  both  conventional  and  incoherent 
(envelope)  beams 

-  A  linear  STA/LTA  detector  applied  to  each  beam 

Frequency-wavenumber  analysis  of  detected  signals 
Association  of  regional  phases  to  aid  in  locating  events. 

In  the  reporting  period  (April-September  1985),  the  RONAPP  system 
has  been  operated  with  a  fixed  beam  deployment,  as  summarized  in 
Table  VIl.3.1.  This  deployment,  which  at  the  current  stage  is  to 
be  considered  experimental,  has  been  based  upon  the  desire  to 
obtain  a  balance  between  teleselsmic  and  regional  phase  detec¬ 
tion,  within  the  limits  of  the  processing  system  capability. 
Briefly,  the  current  deployment  is  described  as  follows: 

Beams  1-7  and  beam  17  are  infinite-velocity  conventional  beams 
with  filter  bands  ranging  from  1-3  Hz  (beam  1)  to  8-16  Hz 
(beam  7).  For  each  of  these  beams,  a  subset  of  the  array  has  been 
selected  as  shown  in  the  table.  Generally,  the  low-frequency 
beams  have  been  formed  from  the  outer  rings  (C  and  U),  whereas 
the  highest  frequency  beams  are  based  on  the  inner  (A  and  B) 
rings.  It  has  previously  been  found  (Mykkeltveit  et  al,  1985) 
that  selecting  such  frequency-dependent  subgeometries  leads  to 


increased  SNR  relative  to  utilizing  all  sensors  of  the  array  in 
the  beamforming. 

Beams  8-16  are  steered  beams  of  typical  teleseismic  velocity 
(14  km/s),  and  variable  azimuths  and  filter  bands.  All  of  these 
beams  are  steered  towards  Eurasia  (azimuth  range  0-180  deg),  with 
the  intention  of  providing  improved  teleseismic  detection  of  the 
high-frequency  signals  seen  from  this  area. 

Beams  18-20  are  Incoherent  (envelope)  beams,  which  have  been 
formed  by  filtering  Individual  sensor  traces  and  then  adding  the 
absolute  values  of  the  filtered  traces  with  no  time  shift 
involved.  As  shown  by  Ringdal  et  al  (1973),  such  beams  are  espe¬ 
cially  well  suited  for  detecting  signals  of  low  coherency  across 
the  array.  This  would,  for  NORESS,  typically  occur  for  high  fre¬ 
quency  secondary  phases. 

We  comment  that  the  RONAPP  package  is  very  flexible  with  respect 
to  beam  deployment,  within  the  overall  limitation  of  20  beams. 
Currently,  the  heavy  filtering  load  involved  in  forming  inco¬ 
herent  beams  limits  the  number  of  such  beams  to  3,  otherwise,  it 
would  have  been  desirable  to  include  more  high-frequency  beams  of 
this  type.  The  threshold  setting  in  the  current  RONAPP  system  has 
been  somewhat  conservative,  in  order  to  avoid  extensive  f-k  pro¬ 
cessing  of  noise  detections.  In  the  near  future,  attempts  will  be 
made  to  allow  for  lower  SNR  thresholds. 

In  the  following,  we  present  some  initial  statistics  on  NORESS 
on-line  detection  results.  These  are  based  on  the  time  period 
April-September  1985,  and  cover  only  those  intervals  (about  75 
per  cent  of  total  time)  when  no  significant  land  line 
transmission  problems  were  encountered.  Such  problems  were 
responsible  both  for  occasional  system  outage  and  periods  of 
large  number  of  spike  detections.  It  is  noted  here  that  the 


satellite  transmission  of  NORESS  data  to  the  U.S.  has  been  con¬ 
ducted  with  minimal  disruption  in  the  period,  thus  the  reliabi¬ 
lity  of  the  NORESS  field  system  has  been  excellent  throughout. 

Fig.  Vll.3.1  shows  NORESS  detections  distributed  by  time  of  day 
(GMT).  The  most  noteworthy  feature  is  the  large  increase  in 
overall  detection  numbers  during  daytime  hours  (local  time  is  1-2 
hours  different  from  GMT).  This  is  clearly  attributable  to  mining 
explosions,  road  work,  etc.,  and  the  figure  shows  that  the  number 
of  detected  local  and  regional  P  phases  (velocity  6-12.5  km/s) 
Increases  in  the  same  proportion  as  the  total  detection  numbers 
during  daytime  hours.  In  contrast,  the  number  of  detections  with 
teleselsmic  P-velocities  (>12.5  km/s)  remains  nearly  independent 
of  time  of  day.  With  the  current  parameter  setting,  there  has  in 
average  been  102  NORESS  detections  per  data-day,  of  which  25  are 
of  estimated  velocity  6-12.5  km/s  and  18  of  velocity  >12.5  km/s. 

In  Fig.  VH.3.2,  NORESS  detections  are  plotted  by  beam  number. 
Note  that  in  this,  as  well  as  the  other  statistics,  only  the  beam 
with  largest  SNR  is  counted  for  each  detection,  and  that  the 
RONAPP  processor  groups  detections  together  if  they  are  separated 
in  time  by  less  than  4  seconds.  From  the  figure  we  see  that  most 
of  the  detections  with  teleselsmic  phase  velocities  occur  on 
beams  1-5,  whereas  the  Incoherent  beams  (18-20)  represent  almost 
exclusively  detections  with  phase  velocities  <6  km/s  (typical 
for  secondary  phases).  A  significant  number  of  detections  with 
low  and  intermediate  phase  velocities  also  occur  on  the  high- 
frequency  infinite  velocity  beams  5-7. 

The  distribution  of  NORESS  detections  by  estimated  azimuth  is 
shown  in  Fig.  VII. 3. 3.  The  concentration  of  detections  in  the 
azimuth  Interval  90-120  degrees  is  mostly  associated  with  mining 
activity  in  central  Sweden  and  western  Russia.  With  regard  to 
detections  of  estimated  teleselsmic  velocity,  the  large  majority 
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is  in  the  range  0-120  degrees,  covering  mostly  Asia.  It  should  be 
noted  that  both  the  azimuth  and  velocity  estimates  computed  by 
RONAPP  have  a  fair  amount  of  uncertainty,  which  must  be  taken 
into  consideration.  This  will  be  further  discussed  in  subsection 
VII. 4. 

Fig.  VII. 3. 4  shows  incremental  and  cumulative  recurrence  sta¬ 
tistics  for  NORESS  detections  as  a  function  of  SNR.  Coherent  and 
incoherent  (envelope)  detections  are  plotted  separately.  For 
coherent  beam  detections,  log  N  (N  =  number  of  detections) 
increases  approximately  linearly  with  decreasing  log(SNR).  On  the 
other  hand,  the  number  of  incoherent  beam  detections  shows  a 
sharp  Increase  at  low  SNR  values.  This  is  partly  because  SNR  of 
secondary  phases  almost  always  is  very  low,  since  the  "noise" 
preceding  such  phases  is  the  P-coda  from  the  same  event.  However, 
preliminary  review  has  also  shown  that  the  number  of  "noise" 
detections  is  significant  for  the  envelope  beams.  Many  of  these 
have  a  very  low  apparent  velocity  (<  3  km/s),  and  procedures  to 
eliminate  such  false  alarms  are  currently  being  investigated. 

In  Fig.  VII. 3. 5,  we  have  subdivided  the  coherent  beam  detections 
into  the  three  velocity  intervals  described  earlier,  and  plotted 
separate  recurrence  statistics.  For  velocities  <  6  km/s,  the 
number  of  detections  increases  very  sharply  with  decreasing  SNR, 
again  probably  due  to  a  predominance  of  low  SNR  secondary  phases. 
(For  coherent  beams,  there  appears  to  be  few  instances  of  "noise" 
detections  due  to  a  conservative  SNR  setting.)  The  plots  for 
intermediate  and  high  velocity  signals  are  similar,  although  the 
"slope"  in  the  latter  case  is  somewhat  less  steep.  In  either 
case,  we  do  not  see  a  sharp  Increase  in  the  number  of  detections 
at  low  SNR,  and  this  is  consistent  with  the  above  assertion  that 
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"noise"  detections  are  essentially  absent  with  the  current  RONAPP 
parameter  setting. 

In  conclusion,  the  six-month  period  of  RONAPP  operation  with 
fixed  parameter  setting  has  provided  valuable  statistical  infor¬ 
mation  which  will  now  be  used  in  incorporating  improvements  into 
the  processing  system.  In  parallel  with  this,  results  are  now 
being  obtained  from  a  systematic  review  of  NORESS  bulletin  output 
which  has  recently  been  initiated.  Altogether,  we  plan,  in  the 
near  future,  to  generate  a  new  beam  deployment.  Including 
possible  lowering  of  the  coherent  beam  SNR  thresholds.  The  intro¬ 
duction  of  variable  thresholds  for  Incoherent  beams  (lower 
thresholds  immediately  following  regional  P  detections)  will  be 
considered  as  a  possibility  to  improve  secondary  phase  detection, 
as  will  the  possibility  of  masking  "noise"  detections  with  very 
slow  phase  velocities.  Another  planned  refinement  is  to  introduce 
a  subroutine  for  faster  f-k  analysis,  thus  eliminating  a  serious 
bottleneck  in  the  real-time  processing,  and  to  implement  an 
interpolation  procedure  to  obtain  better  precision  in  phase  velo¬ 
city  and  azimuth  estimates  than  the  current  grid  points  allow. 
Finally,  the  regional  signal  characteristics  discussed  in  subsec¬ 
tions  VII. 7,1  and  VII. 7. 2  will  be  taken  into  account  in  improving 
automatic  phase  association  procedures. 

F.  Ringdal 
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*********  NORESS  ON-LINE  BEAM  DEPLOYMENT  APR-SEPT  1985  ********* 
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Table  VII. 3.1  NORESS  on-line  processor  beam  deployment  effective 
14  March  1985. 
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Fig.  VII. 3.1  Distribution  of  NORESS  detections  as  a  function  of 
time  of  day  (GMT).  The  average  number  of  detections 
per  hour  is  plotted,  based  on  the  time  interval 

April-September  1985,  with  correction  applied  for  ' 

system  down  times.  The  total  heights  of  the  columns  j 

correspond  to  ail  detections  within  the  respective  _ i 

hours.  ..'j 
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velocity  <  b  km/s. 
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velocity  between  b  and  12.5  km/s,  ;.0 

Doubly  hatched  areas ;  detections  with  phase  velo-  _ 

city  exceeding  12.5  km/s. 
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Fig.  VII. 3. 2  NORESS  detections  during  Aprll-September  1985 

distributed  by  beam  number  (see  Table  VII.3.1).  For 
each  detected  phase,  the  beam  with  the  Largest  SNR 
has  been  selected.  The  hatching  convention  is  the 
same  as  in  Fig.  VII.3.1. 
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NORESS  DETECTIONS  BY  AZIMUTH 
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Fig.  VII. 3. 3  NORESS  detections  during  Aprll-September  1985  as  a 
function  of  estimated  azimuth.  The  hatching  conven¬ 
tion  is  the  same  as  in  Fig.  VII.3.1. 
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Fig.  VII. 3. 4  NORESS  detection  statistics  during  Aprii-September  1985  as  a  function  of  SNR.  Incremental  sta¬ 
tistics  (histograms)  and  cumulative  statistics  (lines)  are  plotted  separately  for  all  detections 
on  conventional  beams  (nos.  1-17)  and  envelope  beams  (18-20)  defined  in  Table  VII. 3.1.  Note  the 
differences  in  horizontal  scales  of  the  two  plots. 
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a  study  has  been  conducted  comparing  the  NORESS  detection  file  to 
the  seismic  bulletin  generated  from  NORSAR  data.  Such  a  com¬ 
parison  is  of  particular  interest  since  the  small  NORESS  array  is 
located  in  the  same  geographical  area  as  the  large  NORSAR  array 
(Fig.  VH.4.1).  Thus  we  obtain  a  measure  of  how  the  high- 
frequency  NORESS  processing  based  on  densely  deployed  sensors 
compares  to  the  NORSAR  processing,  which  utilizes  lower  signal 
frequencies  and  larger  intersensor  spacing. 


The  data  base  for  this  study  covers  the  six-month  period  April- 
September  1985,  during  which  the  NORESS  experimental  detection 
processor  (RONAPP)  has  been  run  with  a  fixed  parameter  setting 
(reference  subsection  VII. 3).  For  the  purposes  of  this  evalu¬ 
ation,  time  Intervals  when  NORESS  operation  was  degraded  due  to 
transmission  line  problems  have  been  deleted,  thus  the  total 
number  of  full  data  days  is  138,  or  approximately  75  per  cent  of 
total  time. 


The  NORSAR  event  list  for  these  138  days  comprises  in  total  1628 
seismic  events.  Ail  of  these  have  been  reviewed  and  accepted  by 
the  analyst,  and  thus  constitutes  an  excellent  reference  data 
base.  The  NORESS  automatic  detection  list  for  the  same  Interval 
contains  14069  entries,  and  has  not  been  subject  to  systematic 
analyst  review.  From  the  results  of  subsection  V11.3,  it  is  clear 
that  the  large  majority  of  NORESS  detections  correspond  to  local 
and  regional  phases,  which  are  not  reported  in  the  NORSAR 
bulletin.  Thus,  the  respective  numbers  of  NORSAR-NORESS  detections 


alone  give  no  measure  of  the  relative  array  detection  perfor¬ 
mance,  and  it  is  necessary  to  match  the  entries  directly  against 
one  another  in  order  to  obtain  a  comparison. 

The  following  procedure  was  applied  in  matching  NORESS  detection 
entries  against  NORSAR-reported  events.  For  each  NORSAR  event, 
the  predicted  P  arrival  time  at  NORESS  was  calculated  based  on 
the  NORSAR  arrival  time  and  the  NORSAR  estimated  phase  velocity 
and  azimuth.  If  one  or  more  NORESS  detections  were  found  within 
15  seconds  of  the  predicted  time,  the  time  of  the  earliest  such 
detection  was  then  examined.  If  this  differed  less  than  3 
seconds  from  the  predicted  time,  a  match  was  automatically 
declared.  If  the  difference  was  in  the  range  3-15  seconds  (about 
5  per  cent  of  the  cases),  the  NORESS  recordings  were  reviewed  in 
order  to  determine  whether  or  not  they  represented  the  same  event 
(typically,  this  would  be  a  emergent  P  phase).  In  cases  where  no 
NORESS  entry  was  found  within  15  seconds  of  predicted  time, 
the  event  was  considered  not  detected  by  NORESS. 

A  total  of  1376  of  the  reference  events  was  thus  found  to  be 
detected  at  NORESS,  i.e.,  8A.5  per  cent.  The  NORSAR-reported 
events  are  almost  exclusively  at  teleseismic  and  greater  distan¬ 
ces,  thus  this  high  percentage  shows  that  NORESS  has  an  excellent 
detection  capability  for  teleseismic  and  core  phases.  It  might  be 
added  that  a  significant  number  of  the  NORESS  detections  which 
have  no  matching  NORSAR  entry  also  appear  to  be  of  teleseismic 
origin,  and  this  will  be  the  subject  of  further  study  in  the 
future. 

Fig.  Vll.4.2  illustrates  the  ditferences  in  observed  and  pre¬ 
dicted  NORESS  arrival  time  for  the  1376  common  events.  In 
reviewing  this  figure,  it  must  be  kept  in  mind  that  the  NORSAR 
arrival  times  have  been  refined  through  analyst  review,  whereas 
the  NORESS  times  have  been  determined  automatically.  There  is  on 


the  average  a  slight  positive  bias  (about  0.5  seconds)  in  the 
data,  and  about  80  per  cent  of  the  observed  differences  are 
within  1.5  seconds  of  this  average.  This  indicates  that  the  auto¬ 
matic  arrival  time  definition  used  in  the  NORESS  processing 
system  works  quite  well,  especially  taking  into  account  that  the 
large  majority  of  the  events  are  of  low  magnitude,  and  many  are 
characterized  by  emerging  P-waves  or  multiple  phases. 

The  percentages  of  NORSAR  events  detected  by  NORESS  are  shown  as 
a  function  of  distance  and  azimuth  (NORSAR  estimates)  in 
Figs.  VII. 4. 3  and  VII. 4. 4.  Not  surprisingly,  NORESS  detects  all 
of  the  NORSAR  events  within  20  degrees  distance  (as  well  as  many 
more  events  not  visible  on  NORSAR  recordings).  Otherwise, 

Fig.  VII. 4. 3  shows  that  the  relative  NORESS  detection  performance 
of  P-waves  decreases  slightly  with  increasing  distance.  A 
noteworthy  feature  is  the  high  percentage  of  core  phases  detected 
by  NORESS;  this  is  primarily  due  to  the  predominantly  high  fre¬ 
quency  signals  seen  from  the  Kermadec-Fi ji  Islands  region  in 
combination  with  favorable  signal-focusing  effects  at  NORESS  for 
this  region. 

From  Fig.  VII. 4. 4  it  is  seen  that  the  large  majority  of  NORSAR 
events  are  in  the  azimuth  range  0-120  degrees.  The  NORESS  detec¬ 
tion  percentages  are  about  constant  when  averaged  over  30  degrees 
azimuth  windows.  The  only  exception,  the  interval  180-210 
degrees,  comprises  very  few  events,  mostly  low-frequency  earth¬ 
quakes  from  the  South  Atlantic  Ridge.  Even  when  taking  the  uneven 
distribution  of  reference  events  into  account,  it  might  be 
surprising  that  NORESS  does  not  show  higher  percentages  for  the 
predominantly  high  frequency  signals  from  Asia  (azimuths  0-120 
degrees)  than  for  the  low  frequency  signals  from  North  America 
(azimuths  270-360  degrees).  The  explanation  lies  in  the  signal 
focusing  effects  across  NORSAR:  For  high  seismicity  areas  such  as 
the  Japan-Kamchatka  arc  and  Hindu  Kush,  NORSAR  benefits  from 


sharp  focusing  effects  at  subarrays  02B  and  03C,  respectively, 
whereas  the  NORESS  site  (U6C)  is  relatively  poor  in  comparison. 
The  fact  that  the  detection  percentages  remain  so  high  is  attri¬ 
buted  to  the  improved  gains  by  the  high  frequency  signal  pro¬ 
cessing  at  NORESS.  The  high  NORESS  detection  percentages  from 
North  America  are  due  to  relatively  favorable  focusing  effects  at 
the  NORESS  site  for  this  region,  which  at  least  partly  compensate 
for  the  detectability  loss  caused  by  low  signal  frequencies. 

For  those  teleseismic  regions  that  are  characterized  by  both 
high-frequency  signals  and  favorable  signal  focusing  effects  at 
the  NORESS  site,  the  NORESS  detection  performance  appears  to  be 
significantly  better  than  that  of  NORSAR  (judging  from  SNR  of 
common  events).  For  example,  Semipalatinsk  explosions  show  10-14 
dB  higher  SNR  for  NORESS  when  comparing  the  respective  array 
beams.  Admittedly,  at  least  part  of  this  difference  will  be  com¬ 
pensated  for  by  the  introduction  of  higher  frequency  filters  at 
NORSAR,  which  is  now  in  progress. 

The  1376  common  NORSAK-NORESS  events  have  also  provided  a  data 
base  for  comparing  various  signal  parameter  estimates. 

Fig.  VII. 4. 5  shows  the  results  with  respect  to  azimuth  and 
slowness.  The  median  NORSAR-NORESS  difference  is  12  degrees 
(azimuth)  and  1.5  s/deg  (slowness).  It  should  be  noted  that  the 
grid  used  in  the  automatic  f-k  analysis  of  NORESS  signals  is 
relatively  coarse,  and  this  is  likely  the  reason  for  some  of  the 
differences.  Also,  at  low  SNR,  the  NORESS  f-k  solution  becomes 
more  unreliable  because  of  noise  interference,  and  analyst 
interaction  would  be  necessary  to  improve  the  estimates  in  such 
cases.  Nevertheless,  it  appears  that  NORESS  can  provide,  in  most 
cases,  a  useful  automatic  estimate  of  teleseismic  phase  velocity 
and  azimuth.  The  eventual  capabilities  of  the  array  in  this 
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regard  cannot  be  assessed  until  more  detailed  regional  studies 
have  been  conducted. 

The  differences  in  estimated  dominant  signal  period  for  NORSAK 
and  NOilESS  are  illustrated  in  Fig.  VII. 4. 6.  This  parameter  is 
measured  manually  at  NORSAR,  automatically  at  NORESS,  in  both 
cases  on  the  respective  array  beam.  The  NORSAR  beam  is 
unfiltered,  whereas  a  variable  prefilter  is  applied  at  NORESS, 
and  this  accounts  for  some  of  the  differences.  However,  the  most 
significant  factor  appears  to  be  the  array  beam  loss  at  high  fre¬ 
quencies  at  NORSAR,  causing  the  NORSAR  period  measurement  to  be 
biased  high  in  many  cases. 

In  conclusion,  this  study  has  shown  that  the  high-frequency 
signal  processing  conducted  at  NORESS  provides  a  teleseismic 
detection  capability  which  is  comparable  to  that  of  the  much 
larger  aperture  NORSAR  array.  NORESS  also  gives  useful  azimuth 
and  phase  velocity  information  for  teleseismic  events,  although 
its  capabilities  in  this  regard  are  limited  by  the  array  aper¬ 
ture.  Further  studies  will  be  directed  toward  continued  improve¬ 
ments  in  NORESS  automatic  processing,  and  also  toward  taking 
advantage  of  the  NORESS  experience  in  improving  high-frequency 
signal  processing  at  NORSAR. 


F.  Ringdal 
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Fig.  1  Map  showing  the  relative  geometries  of  the  NORSAR  and 

NORESS  array.  The  current  NORSAR  SF  configuration 
comprises  7  subarrays  (large  filled  circles),  each  with 
SPZ  seismometers  deployed  over  an  area  10  km  in 
diameter.  NORESS  comprises  25  SPZ  sensors  within  a  3  km 
diameter  aperture. 
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Fig.  VII. A. 2  Histogram  showing  the  distribution  of  differences  be¬ 
tween  signal  onset  times  computed  for  NORESS  and  pre¬ 
dicted  onset  times  based  on  the  NORSAR  bulletin.  The 
data  base  consists  of  seismic  events  detected  by  both 
systems  during  Aprll-September  1985. 


EPICENTRAL  DISTANCE  (DEG) 


Fig.  VII. 4. 3  Top  part  of  figure  shows  the  distribution  by  epicentral 
distance  (20  degrees  increment)  of  seismic  events 
reported  by  NORSAR  April-September  1985,  during  time 
Intervals  when  NORESS  was  operational.  The  bottom  part 
shows  percentages  of  these  events  detected  by  NORESS. 
Note  the  high  NORESS  detection  percentages  at  all 
distances . 
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Hg.  VII. 4. 4  Same  as  Kig.  VII. 4. 3,  but  with  event  distribution  based 
on  NORSAR  azimuth  (30  deg  increments). 
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Fig.  VII. 4. 5  Histograms  showing  differences  (absolute  values)  between 
NORSAR  and  NORESS  azimuth  and  slowness  estimates  for  the 
common  events  in  the  data  base. 
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fig.  VII. 4. 6  Comparison  of  dominant  signal  periods  estimated  by 

NORSAR  and  NURESS  for  common  events  in  the  data  base. 
The  NORSAR  histograms  have  vertical  hatching,  vs 
diagonal  hatching  for  NORESS.  Note  the  significantly 
lower  dominant  periods  estimated  from  NURESS  data. 
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VII .5 _ P-wave  spectra  from  NORESS  recordings 

A  number  of  earlier  studies,  some  of  which  based  on  NORSAK  data, 
have  shown  that  high  frequency  signal  propagation  in  Eurasia  is 
very  efficient,  even  at  teleseismic  distances.  With  the  new 
NORESS  installation  in  Norway,  we  now  have  the  possibility  to 
study  both  the  spectral  contents  and  spatial  coherency  of  tele- 
seismic  phases  at  higher  frequencies  than  before,  since  the 
NORESS  SP  system,  which  is  based  on  40  Hz  sampling  rate,  effec¬ 
tively  records  signals  of  up  to  about  15  Hz  frequency. 

We  have  initiated  a  systematic  study  of  the  spectral  charac¬ 
teristics  of  NORESS-recorded  seismic  phases,  and  in  this  paper  we 
present  some  representative  examples,  together  with  brief  com¬ 
ments  summarizing  the  most  important  features.  Five  events  of  low 
and  intermediate  magnitude  (m^  4. 0-5.0)  have  been  chosen,  at 
least  three  of  these  are  presumed  underground  explosions.  The 
events  are  in  the  epicentral  distance  range  14-49  degrees.  In 
Figs.  VII. 5.1  to  VII. 5. 5  we  show,  for  each  event,  NORESS  single 
sensor  and  beam  time  domain  traces  as  well  as  spectral  plots  of 
signal  and  noise.  The  beam  has  been  steered  to  the  event  epi¬ 
center  using  plane-wave  time  delays,  and  is  based  on  the  subset 
of  NORESS  instruments  comprising  the  center  instrument  and  the  0 
and  D  rings,  i.e.,  17  of  the  25  NORESS  SPZ  channels.  As  discussed 
by  Mykkeltveit  et  ai  (1985),  this  subset  provides  particularly 
efficient  noise  suppression  in  the  band  1. 5-3.0  Hz,  and  is  better 
than  the  full  array  in  this  regard. 

The  noise  spectra  have  been  estimated  using  the  indirect  co- 
variance  method.  We  first  estimate  the  correlation  function  by 
splitting  a  long  data  record  into  many  windows,  calculating  a 
sample  correlation  function  for  each  window,  then  averaging  the 
sample  correlation  functions.  Typically,  we  use  20  windows,  each 
of  which  is  5  seconds  long.  Because  the  earth  noise  has  such  a 
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large  dynamic  range,  we  prewhiten  it  prior  to  estimating  the 
correlation  function  with  a  low-order  prediction-error  filter. 

The  spectrum  is  then  estimated  by  windowing  the  correlation  func¬ 
tion  with  a  3-second  Hamming  window,  then  computing  the  Fourier 
transform.  The  spectral  estimate  obtained  this  way  is  compensated 
then  for  the  effects  of  prewhitening  and  normalized  to  a  1-second 
window  length. 

Signal  spectra  have  been  estimated  using  the  same  technique,  but 
with  4  overlapping  windows,  each  of  7  seconds  length.  Start  times 
of  these  windows  are  3,  2,  1  and  0  seconds  before  signal  onset, 
respectively.  Thus  we  achieve  a  smoothing  of  the  signal  spectra 
while  retaining  compatibility  with  the  noise  spectra. 

The  NORESS  digital  recording  system  employes  gain-ranged  16  bits 
data  words,  as  described  by  Stokes  (1982).  While  the  overall 
dynamic  range  is  as  large  as  120  dB,  the  actual  resolution  at  any 
given  gain  level  is  considerably  lower.  Thus,  quantization 
effects  become  significant  for  those  portions  of  the  spectra 
which  are  70  dB  or  more  below  the  spectral  peak  of  a  given 
signal,  and  this  must  be  kept  in  mind  when  reviewing  the  spectra 
displayed  in  this  subsection.  The  quantization  effects  are  mani¬ 
fested  in  these  plots  as  an  oscillation  in  the  spectral  traces, 
which  is  apparent  on  some  events  at  very  high  frequencies. 

In  the  following,  we  briefly  comment  on  the  main  features  of  the 
five  events  selected; 

Event  1 :  66. ON  40. 8E,  distance  14  deg,  azimuth  55  deg,  m[j=5.0 

The  P  wave  is  very  complex,  with  significant  high  frequency 
energy,  especially  in  the  P  coda.  Much  of  the  high  frequency 
energy  is  lost  on  the  beam,  as  seen  both  from  the  time  domain 
plots  and  the  spectra.  There  is  significant  '^NR  over  the  entire 
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frequency  band  shown,  and  in  fact  Is  appears  likely  that  there 
would  be  significant  energy  also  above  15-20  Hz  for  this  event. 
However,  the  largest  SNR  is  found  in  the  2-8  Hz  band.  This  event 
is  in  the  regional  distance  range,  and  has  large  Sn  and  Lg  phases 
recorded  by  NORESS. 

Event  2 ;  46. 8N,  48. IE,  distance  25  deg,  azimuth  107  deg,  mt,=5.0 

The  P  wave  is  complex  also  for  this  event,  due  to  multiple  onsets 
from  different  upper  mantle  P  phases.  The  signal  has  a  clear 
spectral  density  peak  at  2  Hz,  and  the  spectral  density  drops 
much  more  rapidly  with  increasing  frequency  than  for  Event  1. 
There  appears  to  be  no  significant  SNR  above  12  Hz  frequency, 
even  though  the  quantization  noise  makes  a  definite  conclusion 
difficult.  Beam  signal  loss  is  negligible  up  to  6  Hz  in  this 
case . 

Event  3;  49. 9N  78. IE,  distance  38  deg,  azimuth  75  deg,  mb*4.9 

This  event  is  characterized  by  a  very  clear,  impulsive  P  signal. 

A  sharp  spectral  peak  is  seen  at  2  Hz  for  the  signal,  and  there 
is  significant  SNR  out  to  at  least  15  Hz.  Beam  signal  loss  is 
negligible  up  to  6  Hz.  Note  the  pronounced  noise  suppression  on 
the  beam  around  2  Hz  (associated  with  the  choice  of  subgeometry 
discussed  before);  this  feature,  which  here  is  particularly 
visible,  can  also  be  Identified  on  the  other  figures  by  close 
inspection. 

Event  4 :  Central  Asia  (no  NEIS  location,  distance  ~ 40  deg, 
azimuth  ~80  deg,  mb=4.0 

This  is  a  low  magnitude  event  from  the  same  general  area  as 
Event  3,  and  many  of  the  same  comments  apply.  The  recordings 
illustrate  in  particular  the  large  SNR  gain  on  the  beam,  espe¬ 
cially  around  2  Hz,  and  it  is  noteworthy  that  even  for  an  event 
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at  such  low  magnitude,  a  significant  SNR  can  be  found  on  the  beam 
up  to  around  10  Hz  frequency. 

Event  5;  41. 7N  88. 4E,  distance  49  deg,  azimuth  76  deg,  mb=4.7 

This  event  has  a  clear,  impulsive  P  phase  recorded  at  NORESS.  In 
contrast  to  previous  event,  the  signal  power  around  2  Hz  is  much 
lower  than  the  spectral  maximum,  thus  the  best  SNR  occurs  at  4  Hz 
frequency.  Again,  the  beam  has  significant  SNR  up  to  about  10  Hz. 

In  conclusion,  the  studies  of  these  and  other  events  have  shown 
that  P  phases  from  intermediate  and  low  magnitude  Eurasian  events 
very  often  contain  signal  energy  well  above  the  noise  level  up  to 
at  least  10  Hz  frequency.  However,  for  teleseismic  events,  the 
SNR  on  the  beam  is  usually  best  at  frequencies  below  4  Hz;  thus 
for  teleseismic  event  detection  purposes,  a  set  of  bandpass 
filters  covering  the  range  1-5  Hz  would  be  expected  to  give  the 
best  performance.  At  the  high  frequency  end,  there  is  clearly  a 
potential  for  additional  SNR  gains  through  beamforming  based  on 
array  subgeometries  involving  only  the  inner  rings,  and  this  will 
be  the  topic  of  a  separate  investigation. 

We  emphasize  that  these  results  are  preliminary,  and  we  plan  to 
continue  our  research  on  this  topic. 


T.  Kvsrna 
F.  Ringdai 
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Fig.  VII  5.1  The  top  part  of  the  figure  shows  NORESS  single  sen¬ 
sor  trace  (instrument  AOZ)  and  NORESS  steered  beam 
(employing  the  17  seismometers  AO,  Cl-7,  Dl-9)  for 
Event  1  described  in  the  text.  The  lower  part  of  the 
figure  shows  power  density  spectra  for  (a)  signal  (AOZ), 
(b)  signal  (beam),  (c)  noise  (AOZ)  and  (d)  noise  (beam). 
Tick  marks  on  the  vertical  axis  correspond  to  10  dB 
power  increments.  The  spectra  have  not  been  corrected 
for  system  response. 
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Fig.  VII. 5. 4  Same  as  Fig.  VII. 5.1,  but  for  Event  4  described  in  the 
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^^8*  VH.5.5  Same  as  Fig.  VII. 5.1,  but  for  Event  5  described  in  the 

text.  Due  to  data  problems,  the  spectra  for  this  event 
are  based  on  shorter  time  windows  than  Events  1-4,  and 
the  spectral  estimates  beyond  12  Hz  are  uncertain. 


VII .6 _ The  calculation  of  elllptlclty  corrections 

For  many  purposes,  Including  that  of  evaluating  the  performance 
of  GDSN  stations  in  source  analysis  studies,  it  is  necessary  to 
calculate  synthetics,  or  Green's  functions,  in  a  spherically  sym 
metric  earth  model  (or  a  combination  of  such  models),  but  with 
the  effects  of  asphericlty  of  structure  included  to  first  order, 
the  effects  are  of  primary  significance  to  the  travel  times  of 
waves,  and  the  most  conspicuous  such  effect  is  that  due  to 
ellipticity.  Since  the  ellipticity  correction  is  to  be  applied 
routinely  to  a  variety  of  phases  and  for  many  source-receiver 
combinations,  it  is  important  that  an  effective  computational 
method  is  available.  Dziewonski  and  Gilbert  (1976)  obtained 
expressions  for  the  ellipticity  correction  in  the  form  of  an 
integral  over  epicentral  distance.  On  the  other  hand,  the  spherl 
cal  travel  times  are  usually  evaluated  much  more  efficiently  as 
an  integral  over  radius.  First  order  perturbations  in  these  tra¬ 
vel  times  can  be  evaluated  similarly;  in  fact,  the  ellipticity 
correction  can  be  obtained  as  a  by-product  of  the  spherical  tra¬ 
vel  time  calculation. 

The  one-way  spherical  travel  time  of  a  ray  between  levels  rg  and 
r;  can  be  expressed: 

r,  d(Ti^-p‘^)i 

T  -  J  — - ^ -  (1) 

ro 

where  p  is  the  ray  parameter,  ri  »  r/v,  and  v'  dv/dr  in  the 
spherical  reference  model.  For  a  turning  ray,  rg  is  the  turning 
point  radius.  Following  Dziewonski  and  Gilbert  (1976),  if  the 
asphericlty  is  given  by 


6v(r,0,ii)  -  -v'(r)  5r(r,e,(p) 
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and  letting 

6r(r,9,(t))/r  =  a(r,0,(t)) 
the  perturbation  In  T  can  be  expressed: 


6T  =  /  T)  a -  ±  [a(Ti  -p  )■*  J  (2) 

ro  1  ^  V  ro 

where  In  the  last  term,  +/-  is  to  be  taken  for  down-/upgoing 
waves.  For  a  complete  ray  the  last  term  is  always  cancelled  by 
contributions  from  subsequent  layers,  unless  rQ  or  r^  represents 
a  velocity  discontinuity  or  an  end-point  of  the  ray. 


Numerically,  T  in  eq.  (1)  is  evaluated  by  interpolating  v  between 
ro  and  r^.  For  example,  for  v  =  a(r/rj;)*’! 


T  - - (3) 

1-b  ro 

Similarly,  if  we  choose  (rj-ro)  sufficiently  small  so  that  a  can 
be  averaged: 


2  2  i  , 

6T  =  (i  -  ±  1)  [a(Ti  -p  (4) 

1-b  ro 

To  substitute  for  o  we  need  the  lateral  position  of  the  ray  at  r^ 
and  ro*  This  information  can  be  obtained  from  the  change  in  epi- 
central  distance: 


Further  simplifications  are  possible  if  o  represents  the 
ellipticity.  As  a  check  on  the  accuracy  of  the  calculations  we 
compared  the  results  to  tabulated  results  (Dziewonski  and 
Gilbert,  1976),  and  Fig.  VII. 6.1  shows  an  example  for  surface 
focus  S  waves.  Considering  the  difference  in  reference  model  and 
possibly  in  the  values  for  ellipticity  (we  used  the  values  of 
Bullen,  1975),  the  agreement  is  quite  satisfactory. 

D.J.  Doornbos 
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Seismic  wave  attenuation  characteristics  at  local  and  regional 
distances  in  Norway  and  surrounding  areas  are  known  to  be  very 
dependent  on  the  local  geology,  with  rapid  regional  variations. 
Still,  very  tew  (if  any)  more  extensive  studies  on  this  problem 
have  been  carried  out  in  this  region,  and  we  therefore  now  report 
on  a  work  aimed  at  analyzing  how  the  quality  factor  Q  appears  in 
the  southern  Norway  region. 


Two  different  methods,  the  coda  decay  method  and  the  method  of 
spectral  ratio,  were  tested  out  and  applied  to  local  events 
recorded  within  the  region  studied.  In  all  cases  the  event- 
station  distances  were  less  than  300  km,  and  the  depth  of  the 
events  were  in  the  range  0-35  km.  Because  observed  Q  in  most 
other  studies  is  found  to  depend  strongly  on  frequency  it  is  com¬ 
monly  expressed  in  the  form  Q=Qo*f“  »  where  Qq  is  Q  at  a 
reference  frequency  (usually  1.0  hz),  and  m  is  a  constant. 


The  use  of  the  coda  decay  method  to  determine  Q  was  first  applied 
by  Aki  &  Chouet  (1975)  and  is  in  principle  used  to  investigate 
how  the  coda  energy  from  a  local  earthquake  recording  decays  with 
time.  The  coda  energy  is  thought  to  consist  primarily  of  energy 
bacli-scattered  from  randomly  distributed  inhomogeneities  and  boun¬ 
daries  in  crust,  and  a  medium  which  generates  coda  well  will  then 
be  less  efficient  as  a  wave  transmitter.  Theoretically,  the  coda 
is  assumed  to  be  described  as  a  homogeneous  pool  of  scattered 
energy  covered  in  an  ellipsoid  with  foci  at  the  event's  hypo- 
center  and  the  scation,  and  the  measured  decay  of  this  energy 
with  Lime  will  reflect  1)  within  the  ellipsoid.  For  different 
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scattering  models  it  has  been  shown  that  the  following  simple 
relation  describes  the  decay  of  coda  energy  with  time: 

A(f  ,t)=S(f  )  •t“'^*exp(-iif  t/Q(f  )  ) 

Here  A(f,t)  is  the  present  coda  amplitude,  S(f)  the  coda  source 
factor,  t  the  time  measured  from  origin,  and  f  is  the  frequency. 
The  dispersion  parameter  u  describes  the  type  of  spreading  that 
occur,  and  is  mostly  given  the  value  1.0  which  reflects  scattered 
body-waves.  Linder  the  assumption  that  u  is  a  constant  and  Q  is  a 
constant  in  each  frequency  band,  Q  can  be  found  by  linear  least 
squares  analysis  for  each  fixed  frequency  f  (see  Fig.  V11.7.i). 

Another  way  to  find  Q  directly  from  observations  is  to  use  the 
method  of  spectral  ratio.  This  method  is  perhaps  more  known,  and 
the  way  it  is  applied  in  this  work  is  to  see  how  the  energy  of 
direct  body-waves  decays  between  two  stations  with  approximately 
equal  azimuth  relative  to  the  source.  By  comparing  the  spectral 
contents  of  the  direct  waves  recorded  at  the  two  stations  ,  an 
estimate  of  the  decay  of  energy  as  a  function  of  frequency  is 
obtained  by  linear  least  squares  analysis  (see  Fig.  VII. 7. 2). 

The  estimated  Q  found  using  this  method  will  describe  the  energy 
loss  along  the  refracted  wave-path  of  the  primary  waves  and  will 
therefore  generally  not  describe  the  same  volume  as  Q  found  from 
coda-waves . 

The  data  used  in  this  analysis  are  a  selection  of  19  local  events 
recorded  during  19b2-83  by  the  temporary  digital  network  of 
southern  Norway  (SNN)  operated  by  NOKSAK.  Additionally,  a  selec¬ 
tion  of  14  local  events  recorded  between  January  and  March  1985 
at  the  new  Western  Norway  Network  (WNN),  which  is  operated  by  the 
Seismological  Observatory  in  Bergen,  were  used.  All  stations  were 
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short  period  vertical  instruments,  and  the  events  used  were  rela¬ 
tively  widely  distributed  south  of  62  degrees  north. 

The  coda  decay  method  were  applied  to  recordings  from  both  net¬ 
works,  but  the  spectral  ratio  method  could  only  be  used  for 
recordings  from  the  SNN-network  because  the  stations  in  the  WNN- 
network  were  not  calibrated  at  the  time  this  work  was  carried 
out . 

To  ensure  that  only  high-quality  data  were  used  in  this  analysis, 
only  events  with  acceptable  signal-to-noise  ratio  were  included. 
Tests  showed  that  the  method  of  coda  decay  was  very  sensitive  due 
to  the  noise  level,  but  a  noise  level  less  than  about  20%  of 
the  signal  amplitude  was  not  critical  when  calculating  Q.  This 
observation  was  then  used  as  a  criterion  under  selection  of  data. 
Fig.  VI 1. 7.1  shows  some  examples  of  chosen  coda  signals. 

In  most  studies  on  this  problem  it  is  found  that  the  quality  fac¬ 
tor  Q  increases  with  depth.  This  result  was  also  found  in  this 
work  in  the  way  that  coda  Q  seems  to  Increase  with  increasing 
length  of  the  investigated  time  window  of  data  used.  The  increase 
in  Q  with  window  length  is  tied  to  the  fact  that  the  later  parts 
of  coda  consists  of  energy  back-scattered  from  deeper  parts  of  the 
lithosphere  (with  higher  Q) .  The  observed  coda  Q  did  not  turn  out 
to  differ  significantly  for  subareas  within  the  southern  Norway 
region,  and  the  averaged  Q-values  for  the  entire  region  were 
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Q-value 

Corr.  coeff. 

No  of  obs. 

Window  (secs) 

55*fl.l5 

0.93 

50 

20 

75*f 1.15 

0.94 

46 

30 

120*fl-09 

0.92 

46 

40 

The  time  window  dependence  for  the  spectral  ratio  method  was  not  as 
critical  as  in  the  coda  decay  method.  The  chosen  window  length  was 
5  secs  and  was  based  on  the  experience  that  this  length  gave  the 
most  stable  results.  The  spectral  ratio  method  was  applied  to  P- 
and  S-waves  separately  and  gave  the  following  results: 

Q-value  Waves  No  of  obs.  Qo  m 

12i*f0.89  P  17  121  ±  90  0.89  ±  0.3 

127*fi*08  S  11  127  ±  70  1.08  ±  0.2 

Some  examples  of  the  spectral  ratio  estimation  of  Q  are  shown  in 
Fig.  VII. 7. 2.  It  must  be  emphasized,  however,  that  the  results 
were  rather  scattered,  and  much  more  data  would  be  needed  to 
reduce  the  rather  wide  error  limits. 

The  results  indicate  that  Q  for  P-waves  generally  is  lower  than  Q 
for  S-waves  within  the  region  studied  (see  Fig.  VII. 7. 3).  Another 
interesting  observation  is  that  the  Q-values  for  S-waves  are 
quite  similar  to  those  obtained  for  coda-waves.  This  is  shown 
expllciteiy  for  coda  Q  over  a  40  secs  time  window.  The  fact  that 
Q  for  coda  waves  in  this  window  are  almost  identical  to  spectral 
ratio  Q  values  for  S-waves  probably  means  that  the  coda  primarily 
consists  of  S-waves,  which  also  is  observed  in  other  regions,  and 


that  coda  in  this  window  samples  approximately  the  same  part  of 
crust  as  the  direct  S-waves. 

The  rate  of  frequency  dependence  of  Q  can  be  thought  to  give  an 
indication  of  the  amount  and  size  of  the  scatterers  that  occur.  A 
scattering  model  described  by  Dainty  (1981)  to  explain  the  fre¬ 
quency  dependence  of  Q  showed  that  observed  Q  can  be  divided  into 
two  different  mechanisms,  intrinsic  Q  caused  by  the  medium's 
deviation  from  complete  elasticity,  and  scattering  dependent  Q 
caused  by  inhomogeneities  the  waves  encounter  during  propagation. 
Intrinsic  Q  is  in  most  works  assumed  to  be  relatively  frequency 
independent,  in  which  case  the  frequency  dependence  of  observed  Q 
must  be  mainly  caused  by  scatterers  or  inhomogeneities  in  the 
crust.  The  apparent  energy  loss  due  to  scattering  will  then  be 
most  critical  for  wave  lengths  comparable  with  the  typical  size 
of  these  scatterers,  which  in  turn  means  that  the  attenuation  of 
energy  at  low  frequencies  must  be  stronger  than  expected  from 
intrinsic  attenuation  only.  This  is  clearly  shown  if  we  compare 
the  result  found  in  this  work  with  the  result  from  Hasegawa 
(1985)  for  the  Canadian  Shield  area  (Fig  VII. 7. 4).  These  two 
regions  both  represent  old  shield  areas  and  should  therefore  not 
be  too  dissimilar  geologically.  It  is  seen  from  Fig.  VH.7.4 
that  the  frequency  dependence  of  Q  in  southern  Norway  is  much 
higher,  but  that  Q  for  both  areas  seems  to  coincide  for  higher 
frequencies.  This  can  be  interpreted  in  terms  of  similar  intrin¬ 
sic  0  in  the  two  regions,  but  with  more  and  stronger  inhomoge¬ 
neities  in  southern  Norway. 

Leif  B.  Kvamme 

Jens  Havskov  (Univ.  of  Bergen) 
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Fig  VII. 7.1  Examples  of  coda-recordings  for  one  specific  sta¬ 
tion  band-pass  filtered  at  different  frequencies. 
The  first  part  of  the  recordings  indicate  the  noise 
level  before  P-wave  arrival,  and  on  top  of  the 
signals  the  RMS-envelope  is  dravm  as  an  approxi¬ 
mation  to  A(f,t).  Axes  are  time  after  origin 
versus  amplitude. 
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Fig.  VII, 7. 2  Examples  of  spectral  ratio  estimation  of  Q  for 
P-waves.  On  top  the  Q(f )-f unctions  of  frequency 
obtained  from  spectral  ratio  analysis  of  the  two 
below  recorded  amplitude  spectra.  The  least  squares 
lines  show  the  curve  of  the  best  fit  frequency 
dependent  Q.  The  black  zones  on  the  amplitude 
spectra  Indicate  the  noise  before  P-wave  arrival. 
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Fig.  VII. 7. 4 


Q  as  a  function  of  frequency  for  different  regions 
The  value  of  the  scattering  index  m  is  shown  on  th 
figure.  The  curves  intersect  the  vertical  axis  for 
1  Hz  at  the  value  of  the  proportionality  factor 


The  individual  seismograph  stations  of  global  networks  always 
include  three-component  Instrumentation  for  the  very  simple 
reason  that  the  seismic  wavefieid  comprises  both  vertical  and 
horizontal  ground  motions.  Seismologists  have  for  may  years  suc¬ 
cessfully  exploited  the  information  potential  of  three-component 
records  for  wave  propagation  modelling,  retrieval  of  structural 
Information  (imaging)  ans  source  parameters,  but  these  efforts 
have  mainly  been  confined  to  the  low  frequency  (less  than  0.2  Hz) 
part  of  the  wavefieid.  Much  effort  has  been  invested  in 
extracting  similar  information  in  the  high  frequency  band,  say 
1-10  Hz,  but  efforts  generally  have  met  with  little  success.  The 
reason  for  this  appears  Co  be  twofold: 

1)  High  frequency  records  are  rather  complex  due  to  scat¬ 
tering,  mode  conversions  and  multipathing. 

2)  The  techniques  of  analysis  often  fail  to  produce  wave- 
field  parameters  in  an  easily  interpretable  format.  For 
example,  a  common  procedure  is  to  produce  many  particle- 
motion  plots  which  generally  are  difficult  to  interpret. 

In  this  section  we  present  a  new  approach  for  extracting  wave- 
field  parameters,  based  on  a  priori  models  for  P  and  S  wave  par¬ 
ticle  motions.  Special  attention  has  been  given  to  the  problem  of 
presenting  results  in  an  easily  interpretable  manner  for 
extracting  signal  parameters  convenient  for  a  wide  variety  of 
research  applications.  This  has  provided  some  remarkable  results 
using  data  from  the  NORESS  array,  which  were  presented  in  Che 
previous  Semiannual  Technical  Summary  (Christoff ersson  et  al, 
1985). 


structure  of  3-compoaeat  reglstratloas 

Both  P  and  S  waves  exhibit  a  high  degree  of  linear  polarization. 
Ground  particle  motion  coincides  with  the  azimuth  of  propagation 
for  the  P  phases.  Surface  waves  of  the  Rayleigh  type  are 
generally  elllptlcally  polarized  In  the  vertical-radial  plane, 
the  fundamental  modes  displaying  retrograde  particle  motion  and 
the  higher  modes  prograde  ellipticity.  Surface  Love  waves  are 
also  found  to  be  rectlllnearly  polarized,  but  in  a  horizontal 
plane  orthogonal  to  the  direction  of  wave  propagation.  Micro- 
seismic  background  is  of  the  Rayleigh  type,  but  with  little  pre¬ 
ferred  directionality.  Signal-generated  noise  may  also  be 
polarized,  although  the  direction  of  polarization  is  often  random 
in  nature.  Using  these  various  characteristics  of  polarized  par¬ 
ticle  motion  trajectories,  we  proceed  in  the  following  way: 

Notations ; 

y(b)  =  [yi(t),  y2(t).  y3(t)]* 

is  the  observed  3-component  data  vector  where  yj^  is  the  north 
component,  y2  the  east  and  y3  the  vertical.  The  asterisk  denotes 
the  transpose. 


z(t)  =  [zi(t),  Z2(t),  ...  Z(^(t)]* 


is  a  k-dimensioned  representation  of  the  signal. 
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Xll  ... 


^31  ^3k 

is  a  matrix  of  unknown  constants  relating  the  signal  to  the 
observed  3-component  data. 

e(t)  =  [ei(t),  esCt)]* 

is  the  noise  vector. 

Model 

We  assume  that  the  observed  data  vector  y(t)  allows  the  following 
representation 

yCt)  =  A(t)  e(t)  (1) 

with  the  following  basic  assumptions: 

i)  z(t)  and  e(t)  are  orthogonal  (uncorrelated) 

i:  )  e(t)  has  expectation  zero 

ill)  the  components  of  z(t)  are  linearly  Independent,  i.e., 
the  signal  is  k-dimensional 

tv)  all  moment  up  to  at  least  second  order  exist. 

From  these  assumptions  it  follows  that  the  zero-lag  second  order 
moments  of  y(t)  can  be  written 


E(t)  =  A  (t>(t)  A*  +  'i'(t) 


(2) 


where  E(t)  =  E  [y(t)  y(t)*j 

(!)(t)  =  E  [z(t)  z(t)*] 

'F(t)  =  E  [e(t)  e(t)*] 


Integrating  (or  averaging  in  discrete  time)  eq.  (2)  over  a  time 
window  Tq  <  t  <  T  we  obtain  the  following,  second-order  structure 

Z  =  A  (f  A*  +  '{'  (3) 


E  = 


/ 

To 


Z(t)  dt  = 


1 

T-Tq+I 


T 

1  M 

t=To 


0=  / 
To 


({i(t)  dt  = 


1 

T-Tq+I 


T 

I  4>(t) 

t=To 


T 

/ 

To 


'l'(t)  dt  = 


1 

T-To+1 


T 

I  tU) 

t=To 


Eq.  (3)  is  the  ordinary  factor  analysis  model. 

Before  discusisng  the  estimation  of  the  unknown  parameters, 
usuaily  based  on  the  second  order  measurements,  some  details  on 
the  P  and  S  wave  representations  will  be  given  -  at  this  time  the 
surface  wave  representations  are  excluded. 


Body  wave  representations 

For  P-waves  the  signal  is  1-dimensional  and  will  be  written  as: 
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yi(f-)  =  Zi(t)  +  e^Ct) 
y2(c)  =  \2i  zi(t)  +  e2(t) 

=  ^31  zjCt)  + 

Ln  terms  of  the  second  order  structure  we  have 


=  signal  power  =(1x1  matrix) 


covariance  matrix  of  the  noise. 


For  P-waves  the  second  order  moments  would  be: 


^^o  ^p*  +  'F 


where  (.given  rotated  in  z,  T,  R) 


Ap  =  jo  A  and  't>o  =  (i'll) 

(  \3i  ) 

and  sgn(\ii)  *  sgn(A3i). 


In  the  case  of  S-waves  the  signal  is  2-dimensional,  with  SH  and 
SV  components. 
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yi(t)  =  \ii  zi(t)  +  Z2(t)  +  ei(t) 

y2(t)  =  X.21  z^Ct)  +  X22  Z2(t)  +  e2(t) 

y3(t)  =  0  •  zi(t)  +  \32  Z2(t)  +  e3(t) 


(6) 


where  z^Ct)  is  the  SH  component  and  Z2(t)  the  SV  component. 
In  terms  of  the  second  order  structure  we  have 


^11 

X12 

power  of  SH  0 

A  = 

A2I 

K22 

and 

<))  = 

0 

A32 

L)  power  of  SV 

with  the  columns  of  z  being  orthogonal. 

In  this  case,  the  second  order  moments  would  be: 

Zs  =  As  't’s  As*  +  <?  (7) 

where  (given  rotated  z,  T,  R) 


and  sgn(\i2)  =  sgn(\32). 

The  first  column  in  As  is  the  SH  component  with  signal  power 
())2I,  the  second  is  the  SV  component  with  signal  power  1^22  • 

Estimation 


For  a  3-component  station,  the  E's  in  eqs.  (5)  and  (7)  are  sym¬ 
metric  3x3  matrices  having  at  most  6  different  elements.  In  order 
to  make  the  covariance  structures  identified,  we  can  have  at  most 


6  unknown  parameters  In  A,  and  W*  so  we  must  Introduce  some 
restrictions.  Firstly,  we  will  normalize  the  non-zero  element  in 
the  (t)-matrices  to  be  equal  to  1  which  is  actually  not  a  restric¬ 
tion.  Secondly,  we  assume  that  the  ^-matrix  is  diagonal,  i.e., 
the  noise  is  assumed  to  have  no  zero-lag  correlation.  Further,  in 
some  cases  we  will  use  the  additional  restriction  that  ^  is  pro¬ 
portional  to  some  known  diagonal  matrix  ?o>  “  9  ^0* 

The  estimation  of  the  unknown  parameters  will  be  based  on  the 
observed  second  order  moment  of  the  observations  over  a  time 
window 

1  T 

S  -  77-7  I  y(t)  tCt)*  (8) 

T-Tq+I  t=To 

There  are  several  possible  estimators  that  can  be  used.  Most  are 
based  on  some  fitting  function,  i.e.,  minimizing  the  difference 
between  the  observed  S  and  the  theoretical  E. 

The  most  commonly  used  are 

1)  ML  (MAXIMUM  LIKELIHOOD) 

This  estimator  is  derived  from  Gaussian  assumptions  and 
minimizes : 

F  =  log  |E|  +  tr(SE-l)  -  log  \S\  -  q  (9) 

where  q  is  the  dimension  of  the  data  (in  this  case  q^S), 
and  tr  =  the  trace  operator. 
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il)  ULS  (UNWEIGHTED  LEAST  SQUARES) 

This  estimator  minimizes: 

F  =  trCS-Z)^  (10) 

l.e.,  the  sum  of  squares  of  all  the  elements  In  S-E. 
Ill)  GLS  (GENERALIZED  LEAST  SQUARES) 

This  estimator  minimizes: 

F  =  tr(I-S~lE)2  (11) 

Iv)  PC  (PRINCIPAL  COMPONENT) 

This  estimator  maximizes,  subject  to  normalizing 
conditions  on  A: 

F  =  tr(A*SA)  (12) 


Properties  of  estimators 

The  PC  estimator  leads  to  an  eigenvalue  decomposition  of  the 
observed  moment  matrix  S  and  was  suggested  by  Husebye  et  al  (1967). 
It  differs  from  the  other  three  mainly  In  that  It  focuses  heavily 
on  the  diagonal  of  S  whereas  the  other  put  more  weight  on  the 
of f-dlagonal.  Under  Gaussian's  assumption,  the  ML  and  GLS  estima¬ 
tes  are  asymptotically  equivalent  and  It  Is  possible  to  obtain  a 
statistical  test  for  the  model  fit  to  data.  ML  and  GLS  need  a 
positive  definite  S-matrlx,  whereas  ULS  works  even  with  non- 
grammlan  S.  This  property  would  point  In  favor  of  ULS  because  for 
strong  and  clear  P-waves  much  above  ground  noise,  the  observed  S- 


matrix  will  often  be  near  singular.  However,  this  problem  can  be 
overcome  by  attenuating  the  transverse  component,  l.e.,  by  adding 
white  noise.  We  will  therefore  use  the  ML  estimator,  mainly 
because  of  Its  possibility  for  model  tests.  These  are  of  course 
strictly  valid  only  under  the  Gaussian  assumptions.  However,  as 
has  been  demonstrated  In  the  applications,  the  tests  can  be 
modified  to  work  well  In  practice. 

Model  test 

Under  the  assumption  of  Independent  observations  and  Gaussian 
distribution,  model  tests  can  be  carried  out  based  on  the  ML 
estimator.  We  have  approximately  that 

(1^-1)  •  F  with  =  (T-Tq+I)  (13) 

Is  distributed  as  chi-square  with  degrees  of  freedom  equal  to 
q(q+l)/2  -  number  of  unknown  parameters.  Test  for  the  P-wave  can 
then  be  constructed  In  the  following  way; 

2 

We  assume  that  ¥=((11.  Let  Fj^  and  Fp  denote  the  minimum  of  the 
fitting  functions  for  the  Noise  and  P  models,  respectively,  using 
ML. 

Then  (N^-l)Fjij  Is  approximately  chi-square  with  5  d.f.,  (N^-l)Fp 
approximately  chi-square  with  3  d.f.,  and  the  difference 

(N^~l ) (Ff^Fp)  Is  chi-square  with  2  d.f.  Now,  If  data  really  are 

2 

noise  with  'F  =  <(i  I,  both  the  Noise  and  the  P  model  would  have 
reasonable  fit  to  data  and  the  drop  In  chi-square,  l.e., 

(N^~l ) (Fjj-Fp)  would  not  be  significant. 


For  data  generated  by  a  P-wave,  on  the  other  hand,  the  noise 
model  would  have  poor  fit,  the  drop  in  chi-square  would  be  signi¬ 
ficant  and  the  P  model  would  have  good  fit  to  data. 

Define 

P(P)  =  Probability  x^(3)  >  (N*-l)  Fp 

P(N)  =  "  X^(5)  >  (N*-l)  Ffj  (14) 

P(D)  =  "  X^(2)  >  (N*-l)  (Fp-Fn) 

We  then  use 

P(P)  •  P(PE>0)  =  P(P)(1-P(D))  (15) 

i.e.,  the  "probability  of  a  P-wave  multiplied  with  the  probabi¬ 
lity  of  the  P  energy  being  larger  than  zero". 

This  measure  is  rather  sensitive  to  deviation  from  the  Gaussian 
and  independent  assumptions.  Also,  for  large  N  it  is  sensitive  in 
the  sense  that  small  deviations  from  the  P-wave  model  will  easily 
be  detected.  This  problem  can  be  solved  by  multiplying  with  a 
constant  C(0<C<1),  the  size  of  which  depends  on  the  sampling  rate 
and  frequency.  The  actual  value  of  C  can  be  determined  such  that 
P(P)  •  P(PE>0)  is  small  for  noise  data  and  large  for  data  con¬ 
taining  P  signals.  Similar  tests  can  easily  be  determined  for  S-, 
R-  and  L-waves. 

Concluding  remarks 

The  3-component  analysis  techniques  as  described  above  have  been 
tested  on  a  large  variety  of  seismic  recordings  and  a  few 
examples  here  will  be  mentioned.  For  the  broadband  records  of  the 
recent  large  Mexican  earthquake  (unfiltered)  at  HRV  (Harvard, 
Mass.,  USA)  we  easily  identified  phase  arrivals  like  P,  PP,  PPP, 


PcP,  ScP  and  S.  The  corresponding  azimuth  estimates  deviated  only 
around  ±  2  deg  from  the  true  ones  -  thus,  from  a  single  3- 
component  station  record,  a  good  Initial  epicenter  estimate  Is 
feasible.  Very  many  NURESS  records  have  also  been  analyzed,  and 
for  teleselsmlc  events  adequate  estimates  of  the  slowness-vector 
are  obtainable.  In  the  extreme,  deterministic  scattering  contri¬ 
butions  In  the  P-wave  coda  have  been  Identified.  Our  3-component 
analysis  techniques  also  work  quite  well,  particularly  for  P-wave 
motion,  for  vertical  seismic  profiling  (VSP)  records  where  domi¬ 
nant  signal  frequencies  typically  are  50-100  Hz. 

Another  feature  of  this  technique  is  tht  it  works  quite  well  for 
weak  signals,  which  Is  not  exactly  surprising  as  here  we  exploit 
the  wavefleld  structure  in  contrast  to  stacking  techniques  where 
signal  power  variations  over  the  earth's  surface  are  used. 

One  drawback  with  our  3-comp,  analysis  is  pure  noise  triggering, 
which  is  attributed  to  some  P-wave  energy  in  the  background 
noise;  the  so-called  whispering  mantle  effect.  Therefore,  in  an 
on-line  event  detection  context  3-comp,  triggering  has  to  be 
weighted  by  signal  power  estimates  say  on  the  vertical  component 
before  event  presence  is  declared. 

Finally,  we  consider  that  3-comp,  analysis  should  be  potentially 
very  useful  in  an  "expert  system"  context  (e.g. ,  see  Chen,  1983), 
as  in  this  way  a  systematic  wavefleld  decomposition  is  feasible 


even  at  relatively  low  SNRs,  whereas  conventional  "signal  power 
techniques  like  spectral  and  f-k  analysis  become  inefficient. 


A.  Christoff ersson 
(Uppsala  Univ.) 
E.S.  Husebye 
S.F.  Ingate  (MIT) 
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schemes 


The  high-quality  NORESS  recordings  of  events  at  local  and 
regional  distances  exhibit  many  complex  features  the  significance 
of  which  is  difficult  to  assess  by  visual  inspections  of  the 
records.  What  is  needed  here  is  fast  and  robust  analyzing  schemes 
for  phase-type  identification,  and  estimates  of  the  associated 
slowness  vector.  An  added  advantage  would  be  if  the  above  infor¬ 
mation  or  part  of  it  can  be  used  for  time-domain  filtering  of  the 
original  records  so  as  to  visually  expose  major  features.  The 
ability  to  extract  significantly  more  information  from  local  and 
regional  event  records  is  essential  for  more  advanced  seismogram 
analysis  like  forward  and  Inverse  modelling  and  not  at  least 
expert  system  designs.  In  this  context  we  have  conducted  a  prac¬ 
tical  analysis  experiment  on  NORESS  records  of  local  and  regional 
events,  using  different  techniques,  namely,  f-k,  semblance  and 
the  novel  3-comp,  analysis  technique  described  in  subsection 
VII. 8. 

£-k  analysis 

Our  interest  here  is  to  have  a  computationally  fast,  sliding  win¬ 
dow  approach,  and  the  necessary  details  here  are  as  follows:  Con¬ 
sider  an  array  of  N-sensors  whose  locations  are  rj;  j=l,N.  For  a 
prespecified  time  window  of  the  records,  the  respective  Fourier 
transforms  are  Cj(co).  The  corresponding  frequency-wavenumber 
(f-k)  spectrum  is  defined  as 


_  N  N 

P(oj,k)  =  I  I  Sji(oj)  exp(ik(r j-ri)) 

j*l  1=1 


(1) 


where  Sji  is  the  cross-spectrum  between  the  j-th  and  1-th 
instruments.  Sine  Sji(oj)  =  Cj(w)»Ci(a)),  eq.  (1)  may  be  rewritten 
as 


_  N  N 

P(a),h)  =  I  ),  Cj(u))Ci*((A))  expi(k*rj)  exp-i(k.ri) 

j=l  1=1 


N  _  _  N  _ 

=  I  Cj(w)  expi(k*rj)  •  I  Ci*(a))  •  [expi(k* r^) ] 
j=l  1=1 


=  I  I  Cj(a))  expi(k*rj)|^ 

j=l 


which  only  involves  a  single  summation.  This  technique  for 
P(oj,k)-estimation  was  first  proposed  by  Flinn  and  Smart  (1971). 
Additional  speed-up  in  programming  is  obtained  by  using  proper¬ 
ties  of  the  exponential  function  to  eliminate  indexing  during 
power  component  summation  and  by  using  the  "sum  of  angles"  for¬ 
mula  for  sine  and  cosine  terms  when  estimating  the  Fourier 
spectrum  CjCw).  Frequency  domain  smoothing  was  performed  by  a 
Hanning  operation,  i.e.,  replacing  CjCw)  by 


1  1  1 
—  Cj((i)i)  +  —  Cj(u))  +  —  Cj(a)2) 


(1)]^,  0)2  being  o)  **5  Hz)  during  the  initial  spectrum  estimation. 


Comments :  f-k,  analysis,  including  its  ML  variants,  is  generally 
not  rated  a  robust  analyzing  technique  for  several  reasons.  In  ML 
estimation  the  signal  is  presumed  non-stochastic,  which  is  not 
necessarily  the  case  in  shorter  time  windows,  while  for  conven¬ 
tional  f-k  relatively  long  time  windows  are  needed  (T  >  1.6  sec) 
to  ensure  a  good  Fourier  spectrum  estimate.  Perhaps  the  most 
intriguing  aspect  of  f-k  analysis  is  that  in  most  applications, 
the  slowness  estimate  extracted  is  tied  to  peak  power  in  the  f-k 
space.  Under  adverse  conditions  with  interfering  wavelets,  peak 
power  is  not  always  associated  with  the  "primary"  phase  and  thus 
wild  slowness  estimates  may  ensure.  Let  us  add  that  under  more 
normal  conditions,  conventional  f-k  analysis  represents  in  many 
respects  a  convenient  tool  for  seismogram  decomposition. 

Semblance 

So-called  semblance  analysis  is  in  many  respects  similar  to 
VESPAGRAM  analysis,  and  as  such  is  of  potential  interest  in  pro¬ 
cessing  of  regional  NORESS  records.  The  semblance  parameter  S^, 
essentially  the  normalized  beam  power  over  a  given  window  2AT  and 
for  a  given  slowness,  is  defined  as: 


I  x(p)i 
=  1 


The  St  parameter  is  a  simple  signal  coherency  measure  which 
appears  to  provide  reasonable  phase  velocity  estimates  even  for 
complex  local  and  regional  event  records.  Besides  its  calcula- 
tlonal  simplicity,  an  added  advantage  is  that  the  semblance  tech 
nlque  works  well  even  for  short  time  windows  of  the  order  of  .5 


sec.  In  practical  use  the  looping  over  the  P  and  S  crustal 
velocity  window,  say  3-10  kms~^  for  a  given  azimuth  and  signal 
window  width.  We  may  also  loop  over  azimuth  for  obtaining  a  more 
refined  velocity-azimuth  estimate  for  dominant  phases  in  the 
record. 

3-comp,  signal  analysis 

The  theoretical  basis  for  this  technique  is  described  in  subsec¬ 
tion  VII. 8,  so  here  we  will  only  comment  on  the  practical  usage 
of  this  technique.  The  output  from  the  3-comp,  analysis  (single 
site)  is  probability  of  P,  S-,  Rayleigh  and  Love-waves  in  the 
records  as  a  function  of  time  and  azimuth.  Since  this  technique 
embodies  the  estimation  of  the  principal  axis  of  the  particle 
motion  ellipsoid,  it  is  easy  to  estimate  also  the  apparent  angle 
of  incidence  for  an  incoming  phase.  This  parameter  can  in  turn  be 
converted  to  apparent  velocity  and  henceforth  epicentral  distance 
via  standard  travel  time  tables.  The  error  in  the  corresponding 
azimuth  and  velocity  estimates  are  of  the  order  of  ±  15  deg  and 
±  3  kms~^.  Important,  most  of  the  observed  deviations  are  deter¬ 
ministic,  that  is,  due  to  structural  effects  and  thus  to  a  large 
extent  can  be  removed  by  appropriate  corrections  in  the  same  way 
as  done  for  arrays. 

The  advantages  with  the  novel  3-comp,  analysis  is  that  it  pro¬ 
vides  an  efficient  means  for  decomposition  of  even  complex 
wavetrains,  can  operate  on  short  time  windows  of  the  order  of 
0. 5-1.0  sec  and  provides  good  estimates  of  the  slowness  vector. 

A  drawback  is  easy  triggering  on  noise  wavelets  which  we  consider 
removable  by  appropriate  filtering  schemes. 
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Practical  examples : 

The  relative  merits  of  the  3  analysis  methods  discussed  above 
have  been  tested  on  real  data  and  detailed  results  for  a  presumed 
explosion  near  Leningrad  (Date:  29  Jan  85;  OT:  11.59.47;  Lat.  & 
Long.:  59. 3N,  28. IE)  are  presented  below  in  Figs.  VII. 9. 1-6  and 
Tables  VII. 9. 1-3. 

All  the  figures  have  the  same  time  scale  and  thus  are  easily 
comparable.  The  main  features  are  summarized  in  the  figure 
captions. 

In  Tables  VII. 9. 1-3  the  outputs  from  f-k.,  semblance  and  3-comp, 
analysis  are  presented  in  detail.  All  three  methods  give  reason¬ 
able  velocity  and  azimuth  estimates  of  the  first  part  of  the  P- 
wave  train  (centered  around  t  =  126  sec),  although  the  semblance 
results  by  far  appear  to  be  the  most  robust  ones.  Interestingly, 
the  single  3-comp,  station  analysis  results  contain  all  prin¬ 
cipal  phase  information  features  as  found  in  the  semblance  and 
f-k  results  plus  some  more,  most  of  which  is  termed  coda  signa¬ 
ture  of  the  event  in  question.  Also,  the  3-comp,  analysis  results 
are  rather  stable  when  based  on  beam  traces  formed  from  the  4 
NORESS  stations  in  question.  For  the  shear  wave  section  we  have 
not  run  the  Love  and  Rayleigh  models  of  our  3-comp,  analysis 
package . 

Comments :  The  results  presented  are  rather  self-evident,  so  we 

will  only  remark  that  the  advantage  of  f-k  is  scanning  of  the 
whole  "space",  but  that  in  more  extreme  cases  the  associated  (peak 
power)  slowness  estimate  could  be  very  wrong.  Semblance  is  robust 
and  provides  adequate  slowness  estimates,  but  for  easy  use  in  an 
automated  mode  require  a  priori  information  on  azimuth.  Further¬ 
more,  semblance  function  values  are  convenient  for  generating 
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Output  of  "sliding"  window  f-k  analysis  for  signal  "center" 
frequency  fo*3.2  Hz  and  whole  NORESS  array  except  A-ring. 
Three  different  window  lengths  are  used,  namely:  1.0,  1.5 
and  2.0  sec.  Other  values  of  fQ  gave  less  consistent  results 
vis-3-vls  "true"  velocity  and  azimuth  at  8.38  kms'l  and 
92.60  deg.  A  comparison  with  semblance  results  in  Table 
VII. 9. 2  clearly  Implies  why  semblance  is  termed  a  robust 
analysis  technique  re  f-k. 


FK-*N»LrSI;,  nlli  SLIOIiG  I 


04T»  FILt  ;  'lKSd‘'iO?9 


3.^0 


V-Lu  PuWtR 


2  .OP 
aZ  IM 


ANALYStP  ^,Il'Ii.Y 


rtl  fiPoW  j  I 

1  .30 

T  I  (1^ 

R'"  rtt  K 

AZIM 

2  1  0. OP 

3o.  1 

99.3  9 

2l  0.  3 ' 

0  7.7 

79.3  0 

2  1  i  .  u  '' 

?  tl.  3 

.■'23  .00 

1.  -j: 

9  1.0 

1 01 .32 

2it.u: 

.^6.7 

311.70 

^12. OP 

9  J  .  1 

67.7  1 

21  3.  J-' 

3  0  •  0 

67.. '>.0 

21  3. 3:. 

3  7.9 

2  1  .PZ 

2  1  ^ .  vj 

9..  1 

1  1 '  9 

21  5^ 

09.  1 

1  3P .9  3 

21  3  .  U  J 

99.  1 

1  o7 .9  5 

21  3. 51 

17.3 

133.31 

2  1  a  .  0  0 

J  o.  a 

2  39.09 

21 6. 31 

91.0 

1 33 .Ou 

21  7  .  uO 

9  0.2 

9'+.'1.t 

<-1  7.  30 

93.7 

71  .76 

2 1  d .  0  r 

3  9.9 

3o;>  .2  5 

2 1  a .  30 

9...  1 

112.29 

219. o: 

9o.  3 

121.71 

21 9. 3 ; 

9  2.3 

123.33 

2  2  0.  j 

••3.7 

111.97 

220.30 

9  0  .  9 

66  .99 

22  1.0"' 

9  H  .  0 

o.l  .95 

221. 30 

5  u  .  3 

12.1.1  7 

222.01 

9  0 .  9 

122.9  7 

222.  3"' 

99.  y 

79.39 

223.0' 

9  0 .  3 

393.1  6 

223. 31 

30.9 

7  7  .,62 

22  H. uO 

9  0.2 

111.37 

22  9.  3.' 

9  9.3 

30.99 

225.03 

-7.3 

132  .19 

225.  5") 

9  7.2 

100.31 

226.00 

93.3 

183.33 

22o. 50 

9  w  .  9 

8  6.76 

22  7.  00 

97.9 

99 .09 

227. 50 

92.0 

o3  .92 

2?b.00 

93.7 

69,73 

22  6.  3 

9  :i .  0 

o7.1  1 

229.0^ 

9  2.7 

2  6  1  .  2  7 

229. 30 

9  0  .  (3 

93.3  3 

.,3  0.  uO 

9  9.7 

169.99 

2  3  0.  59 

90.6 

72  .96 

231.0-1 

9  3.0 

123.91 

2  3  i.  3  0 

9  1  . 

1 33  .92 

..32. 01 

Ah.  3 

91 .1  7 

2  3  2.5'' 

9e.5 

lOO .39 

233.00 

99.2 

101.3  1 

233. 3) 

-6.  3 

93.69 

239. oO 

9  7.3 

115.9  7 

239.  31 

99.  0 

9l  .9  7 

233. 00 

9  0 . 0 

69.6  1 

2  3  3.  3  0 

9  3 . 9 

a  .3a 

2  3  t> .  0  0 

9  9.  0 

91  .22 

23o. 3' 

33.9 

91.12 

23  7.  0  3 

‘•0.6 

12.3.3  6 

.’7.3 

'•0.1 

1 93."t 

23o. 00 

*+0.7 

121-30 

23a. 31 

•'  1 .  f 

1  O’ ,9  1 

23  9. w  . 

9o.  2 

77.13 

239. 30 

‘'0.6 

91.17 

29  0. 0  ' 

9  0.  t 

26.30 

3.09 

91.6 

0.  91 

9.3o 

41.2 

1  3 1. 13 

3.9  3 

93.0 

1  1  1. 57 

3.33 

•♦6  .  ! 

105.02 

3.30 

90.9 

79.  51 

9.71 

96  .  3 

11 5.02 

.3.99 

95.. 0 

96.69 

12.0.^ 

93.0 

10^.99 

'..9  7 

92.3 

j3  7. 63 

f  .“3 

91.7 

137.60 

6.62 

93.7 

119. 13 

4.32 

92.0 

1  1  3.  73 

13. ‘'9 

92.2 

91.90 

!<'.  1 1 

9O.8 

95.0O 

9.10 

93.  } 

106.11 

3.37 

99.0 

106.93 

1  3  .  Co 

99.9 

106. 56 

9.97 

96.2 

93.37 

3.89 

<♦5.0 

117.90 

2.67 

99.8 

-  85.81 

3.95 

95.6 

128.98 

12.77 

98.3 

121.37 

7.3  7 

50.8 

83.  66 

9.69 

50.9 

115.19 

6.21 

52.9 

97.82 

9.9© 

53.9 

77.01 

3 , 4i 

99. 0 

393.99 

3.79 

y.i.n 

67.  98 

3.26 

97.9 

120.96 

6.13 

95.0 

40.29 

2.06 

50.1 

96.63 

9.69 

50.6 

68.81 

3.23 

90.7 

126.51 

9.36 

51.3 

8  1.33 

6.33 

30.9 

9b.  13 

4.72 

96.2 

269.29 

9.15 

32.3 

7  3.  39 

9.97 

96.9 

6  1.99 

1.O0 

97.3 

189.  32 

3.16 

53.3 

'12.91 

19.0fJ 

95.2 

1 65. 97 

9.79 

97.9 

7  6.  tt3 

9.32 

98.3 

97. 82 

2.77 

99.9 

99.61 

9.96 

51.9 

88.73 

3.3a 

99.3 

' 118.20 

5.30 

51.3 

98.89 

3.39 

99  •  7 

8.13 

3,60 

96.  7 

279.90 

6.23 

97.2 

258.  18 

1.36 

59.7 

13.99 

9. 9  I 

99.  7 

117.35 

3.  I  7 

51.4 

119.68 

9.76 

53.3 

llY.05 

7.2  3 

53.6 

10  7.  70 

9.  I  i 

9  8.6 

199.09 

9.9o 

99. C 

216. 12 

9  .  ."  6 

5,1.3 

16.01 

9  .  I  6 

50.2 

j2  b.  77 

9 .96 

50. 8 

Zl.  17 

6  .  I  0 

95.9 

2  8.  50 

2.5u 

V5L0 

pn«eK 

AZIM 

Veto 

3.  d6 

9i  .  0 

3.29 

9.5  8 

2 .  j2 

9  6,6 

65.52 

10.56 

5-26 

97.7 

1 07 -  9  5 

6 . 52 

5.72 

47.4 

95.53 

5.5  1 

8.85 

96.2 

112.91 

9.79 

9.69 

49.5 

112.39 

9.97 

9.65 

90.9 

I0I.59 

7.39 

6  .  U7 

9  3.0 

105.26 

’.lo 

6.33 

92 . 9 

335.25 

8.10 

7.80 

9  5.0 

1  8  .  50 

7  .  i1  „ 

7.65 

99  ,  u 

239.16 

:,97 

8.90 

47.  7 

95.  39 

3.  <=6 

5.  93 

9  5 ,  ♦ 

68.15 

7.^16 

9.79 

9  9.  u 

113.70 

9.  =6 

5.  19 

51.3 

I'OI.  3  9 

9.36 

76  .  o9 

50.9 

1  2  1 . 8  0 

3.17 

6.30 

46.8 

2  6  0 .  1  0 

3.?  6 

19.27 

97.2 

90.8  3 

3.52 

9.21 

49.9 

129,56 

9.95 

3.91 

59.  1 

121 .29 

9. 17 

9.88 

52.0 

125.12 

5.3; 

5.58 

99.9 

3.5  2 

3.07 

5.37 

50.6 

75.76 

3.5c 

9.99 

56. ♦ 

107.02 

4.79 

9.72 

52.3 

120.53 

5.37 

6.07 

53.9 

112.01 

9.79 

6.30 

50.8 

127.23 

2.5a 

5.97 

53.9 

71.57 

5.12 

41.67 

51.2 

393.98 

3.16 

9.76 

51.0 

88,73 

5.40 

5.61 

53.0 

I0O.I8 

6.13 

♦  .62 

52.3 

93.37 

4.76 

2.59 

99.9 

128.15 

2.5'y 

♦  .  07 

55.0 

75.22 

9.9j 

6.87 

52.9 

99.06 

3  ,  <*  a 

3.00 

52.5 

267.88 

.3.00 

9.09 

53.9 

71.57 

9.09 

♦  .56 

96.9 

132.10 

12.68 

3.  58 

99.5 

163.81 

3.38 

9.11 

51.6 

63  .  Iq 

17.9*^ 

8.71 

97.2 

1 IG.Ol 

9.56 

S.03 

13.7 

93.37 

9.76 

9.72 

98.7 

3«.0.0o 

3.  3i 

5 . 0  9 

52.2 

96.63 

5.61 

5.90 

5<;  ,  6 

86.76 

9.5a 

3.10 

52.0 

117.53 

3  .C-+ 

5.39 

98.5 

209.55 

3.25 

3h.37 

97.5 

115.18 

3.76 

3.60 

59.6 

91.12 

9.76 

5.53 

5  6.0 

93.37 

9 . 7  0 

5.o9 

57.1 

107.59 

5.5-j 

7.94 

51.  i 

118.56 

5.16 

5.97 

5o  .  6 

120.56 

5.37 

9.72 

59.9 

101.73 

9.99 

9.93 

56.8 

1  0  7 , 7  0 

9.17 

3.66 

53.2 

1h9.09 

3.  .5  6 

5.31 

5  J  . 

2'3.7^ 

1.97 

5.19 

5^.9 

329.97 

5.17 

7.01 

50.2 

326. 77 

7.31 

9.96 

50.8 

91.17 

9.96 

6.10 

45.9 

28.50 

1.10 

Table  VII. 9. lb  Output  of  "sliding"  window  f-k.  analysis  for  shear 
wave  section,  caption  otherwise  as  for  a). 
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Table  VII. 9. 2a  Output  of  "sliding"  window  semblance  analysis  with  azimuth 
fixed  at  90  deg,  and  traces  3. 0-5.0  Hz  bandpass  filtered. 
Windows  of  lengths  0.3,  1.0  and  1.3  have  been  used. 
Semblance  values  less  than  0.3  not  considered  significant. 
Finally,  fixing  azimuths  at  80  and  100  deg  did  not  pro¬ 
foundly  change  the  semblance  results  displayed  here. 
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0.  3b 
0.06 


.00 

YELO 

7.50 

7.50 

7.50 

7.25 

3.75 
4.00 

4.25 

5.25 

5.50 

6.50 

6.75 

6.75 

6.75 

4.25 

4.25 

4.25 
b.OO 
5.00 
7.00 
5.00 

7.50 

5.75 

5.75 
5.  75 

4.75 
4.00 
4.25 
4.  50 

4.50 
4.  50 

4.50 

4.25 

4.25 
5.00 
5.00 
5.00 

5.25 

3.50 

3.75 

3.75 

3.75 

3.75 

7.50 

4.50 
4.  50 

4.50 

4.50 

4.50 
4.  75 
4.  75 

4.75 
4. 7‘ 
4.  50 

4.50 
4.  50 


Table  V11.9.2b  Output  of  "sliding**  window  semblance  analysis  for 
shear  wave  section.  Caption  otherwise  as  for  a). 
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Table  VH.9.3a  Estimates  of  apparent  P-wave  velocity  at  site  Aq  on  the 
basis  of  the  3-comp,  seismogram  analysis.  Window  length 
1.0  sec  and  traces  3. 0-3.0  Hz  bandpass  filtered.  For  P-wave 
presences,  probabilities  less  than  30  per  cent  (0.3),  velo¬ 
city  estimates  deleted  due  to  lack  of  significance. 
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Table  VII. 9. 3b  Estimate  of  apparent  shear  wave  analysis  for  shear 
wave  section;  caption  otherwise  as  for  a). 


Fig.  VII.9.1:  Filtered  3-comp,  traces  for  station  AO;  the  z- 
component  is  representative  for  the  other  NORESS  stations, 

a)  P-wave  section,  filter  3. 0-5.0  Hz;  b)  shear  wave  section, 
filter  2. 0-4.0  Hz. 

Fig.  VII. 9. 2a):  P-wave  semblance  estimates  as  a  function  of  time 
and  phase  velocity  for  T  =  1.0  sec.;  azimuth  90°.  b)  Same  as  a) 
except  that  looping  is  over  azimuth  with  vel  =  9.0  kms“l.  For 
details  see  Table  VII. 9. 2a. 

Fig.  VII. 9. 3a):  Shear  wave  semblance  estimates  as  a  function  of 
time  and  phase  velocity  for  T  =  1.0  sec;  azimuth  90°.  b)  Same  as 

a)  except  that  looping  is  over  azimuth  with  vel  4.5  kms“^.  For 
details  see  Table  VII. 9. 2b. 

Fig.  VII. 9. 4a):  3-component  analysis,  P-wave  probabilities  (0.9  = 
90%)  as  a  function  of  time  and  azimuth;  window  T  =  1.0  sec, 

b)  Same  as  for  a)  except  this  time  for  the  shear  wave  section. 

For  estimates  of  apparent  velocities,  see  Table  VII. 9. 3. 

Fig.  VII, 9, 5a:  Probability  filtered  records;  only  that  part  of  the 
seisaogram  having  probabilities  of  being  P-wave  motion  within  the 
azimuth  sector  90  t  30°  is  retained,  b)  Probability  filtered 
shear  wave  section. 

Fig.  VII.9.6a-e):  This  sequence  of  figures  demonstrates  the 
viability  of  3-corap.  signal  analysis  under  adverse  conditions, 
l.e.,  SNR  around  1  or  less,  a)  Filtered  P-wave  section 
(3. 0-5.0  Hz)  for  the  beam  formed  on  the  basis  of  the  4  NORESS 
3-comp,  stations;  b)  The  probability  filtered  variant  of  a); 

c)  Unfiltered  beam  record;  d)  The  probability  filtered  variant  of 
c)  (the  original  non-band  pass  record);  e)  Band  pass  filtering  of 
the  d)-records,  which  should  be  compared  to  b).  Our  remark  is 
that  this  record  sequence  demonstrates  convincingly  that  the 
3-comp,  analyzing  technique  works  well  for  SNR  around  1  or  even 
smaller  values.  It  should  be  added  that  we  were  unable  to  repro¬ 
duce  a  counterpart  to  e)  on  the  basis  of  a  single  station  3-comp, 
record. 
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The  purpose  of  this  study  is  an  on-going  effort  to  determine  the 
nature  of  the  teleseismic  P  coda  seen  on  short  period  seismo¬ 
grams,  and  specifically  to  separate  it  into  a  contribution  from 
scattering  near  the  source,  if  the  source  is  in  the  crust,  and 
scattering  near  the  receiver.  To  accomplish  this,  digital  data 
from  NORSAR  and  NORESS  were  used  in  the  frequency  range  1-7  Hz 
and  covering  a  time  span  of  20-200  sec  after  first  P.  The 
variation  of  the  coda  power  spectrum  of  various  types  of  events 
with  time  was  examined  and  found  to  fit  a  simple  exponential 
decay  for  all  events,  although  the  amplitude  of  the  coda  relative 
to  first  P  was  quite  different  for  different  types  of  events, 
being  large  for  crustal  earthquakes  and  lower  for  deep-focus 
earthquakes  and  explosions.  Data  from  NORESS  were  examined  using 
frequency-wavenumber  methods,  particularly  by  determining  the 
power  as  a  function  of  wavenumber  at  a  fixed  frequency;  this  is 
equivalent  to  finding  the  power  as  a  function  of  apparent 
velocity.  Deep-focus  events  had  a  coda  dominated  by  low  apparent 
velocities  (3. 5-4. 5  km/sec),  explosions  a  coda  with  equal  power 
in  high  (10  km/sec)  and  low  apparent  velocities,  crustal  earth¬ 
quakes  a  coda  dominated  by  high  apparent  velocities.  These 
results  indicate  that  the  teleseismic  P  coda  in  the  time  and  fre¬ 
quency  range  considered  indeed  consists  of  energy  scattered  near 
the  source  if  the  source  is  in  the  crust  (the  high  apparent  velo¬ 
city  component),  and  energy  scattered  near  the  receiver  (low 
apparent  velocity  component),  and  energy  scattered  near  the 
receiver  (low  apparent  velocity  component).  The  low  velocities 
indicate  that  teleseismic  P  to  Lg  (trapped  shear  wave)  is  the 
dominant  mode  of  scattering  near  the  receiver;  the  differences 
between  crustal  earthquakes  and  explosions  suggest  Lg  to  tele¬ 
seismic  P  near  the  source.  Multiple  scattering  is  probably  oc¬ 
curring  in  the  coda.  These  results  Indicate  that  coda  magnitudes 
will  depend  on  conditions  in  broad  (~  500  km)  regions  around  the 


source  and  receiver,  perhaps  explaining  their  stability.  It  may 
be  possible  to  separate  the  coda  from  the  source  and  receiver 
regions,  leading  to  even  greater  stability.  The  difference  be¬ 
tween  crustal  earthquakes  and  explosions  may  be  useful  as  a 
discriminant.  These  results  are  detailed  in  a  recent  report  by 
Dainty  (1985). 

An  interesting  aspect  of  coda  analysis  is  to  what  extent  scat¬ 
ters  can  be  separated  into  deterministic  and  random  contribu¬ 
tions.  A  first  experiment  here  was  to  identify  clear  crustal 
scatters  (Pg  vel.  <  7  kms*"^,  7.0<Pn<10.0  kms“^)  by  3-comp,  coda 
analysis,  and  then  using  differential  travel  times  (re  P-onset 
time)  to  actually  locate  these  scatters.  Preliminary  results  are 
displayed  in  Figs.  VII. 10.1  and  VII. 10. 2  and  the  following  com¬ 
ments  apply.  Pg-scattering  sources  are  mainly  concentrated  to  the 
north  and  east  of  NORESS,  and  the  apparent  clustering  roughly 
coincides  with  major  tectonic  boundaries.  Surprisingly,  for 
events  to  the  south  hardly  any  crustal  scatters  were  found 
(mainly  non-shield  areas).  Roughly  the  same  results  are  found  for 
Pn-scatters,  but  the  source  locations  here  are  more  dispersed. 

The  scatter  locations  in  this  case  would  be  relatively  less 
constrained,  as  the  scattering  mechanism  may  well  be  a  "bumpy" 
Moho  or  asymmetric  reflections  from  the  surface  (Pn-mode  of  pro¬ 
pagation  would  be  preferential  beyond  say  250  km). 

A.  Dainty  (MIT) 

E.S.  Husebye 
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Fig.  VII. 10.1  Locations  of  secondary  "sources"  for  coda  phases 

of  Pg  wave  type  being  subject  to  scattering  at  the 
receiver  side.  The  thin  lines  represent  azimuth  of 
the  events  used  in  analysis.  Note  that  the  scat¬ 
tering  efficiency  appears  to  be  relatively  weak  to 
the  south  and  west. 
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APPENDIX  A 


NORESS  DATA  TAPE  FORMAT 


The  following  tape  format  has  been  agreed  between  SANDIA,  LIVERMORE 
and  NORSAR*  On  request,  NORSAR  will  provide  NORESS  data  tapes  with  the 
following  options: 


-  1  sec.  blocking  (35A0  bytes  per  record) 

-  2  sec.  blocking  (70b0  bytes  per  record) 

-  10  sec.  blocking  (35400  bytes  per  record)  (standard) 

-  800/1600  or  6250  bpi,  9  track  tapes  (6250  bpi 

is  standard 

-  Array  Geometry  Record  in  EBCDIC  (standard)  or 

ASCII 

A  NORESS  tape  consists  of  3  sections: 

1.  FILE  HEADER  (TAPE  HEADER)  80  bytes 
1  *  EOF 

2.  ARRAY  GEOMETRY  (AG)  RECORD  3600  bytes 
IRG  (INTER  RECORD  GAP) 

3.  DATA  (35  FRAMES  *  100  BYTES  +  RECORD  HEADER  40  BYTES) 

TOTALLY  3540  BYTES  IN  EACH  1  SEC.  BLOCK. 

IRG  between  each  sec. 

Last  1  sec.  block  is  followed  by  EOF. 

Note  that  the  data  record  lengths  will  be  different  if  2  or  10  seconds 
blocking  has  been  requested. 

In  some  cases,  it  might  be  practical  to  include  several  data  intervals 
on  one  magnetic  tape.  These  intervals  will  then  each  comprise  three 
sections  (File  header,  AG  record.  Data),  separated  by  EOF,  unless 
otherwise  specified  in  documentation  accompanying  the  tape. 


1.  FILE  HEADER; 


Bytes  1-14:  Fixed  format  as  follows: 


Bytes 

1-2  : 

2 

byte 

Bytes 

3-4  : 

2 

byte 

Bytes 

5-6  : 

2 

byte 

Bytes 

7-8  : 

2 

byte 

Bytes 

9-10: 

2 

byte 

Bytes 

11-12: 

2 

byte 

Bytes 

13-14: 

2 

byte 

binary ... . Type  00 
01 

binary. . . .Month 
binary. . . .Day 
binary. . . .Year 
binary. . . .Hour 
binary. . . .Min. 
binary. . . . Sec . 


-  AG  record  in  EBCDIC 

-  AG  record  in  ASCII 


This  time  corresponds 
to  the  first  data  on 
the  file 


rv'-  •  -  ' 
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Bytes  15-80:  Variable  format.  Tapes  generated  at  NORSAR 
will  contain: 


byte 

15,16,17,18,19,20 

EBCDIC 

tape  number 

byte 

21,22 

EBCDIC 

Creation  month 

byte 

23,24 

EBCDIC 

Creation  day 

byte 

25-28 

EBCDIC 

Creation  year 

byte 

29-32 

EBCDIC 

Array  Geometry  version  no 

byte 

38-80 

EBCDIC 

comments . 

2.  ARRAY  GEOMETRY:  (EBCDIC  or  ASCII) 
2.1  General  layout 


byte 

1  - 

80 

Text  (free  field)  specifying  date 
and  time  for  last  change  in  array 
geometry.  See  Table  A-1,  first  line 

byte 

81  - 

160 

Text  (Table  A-1,  second  line) 

byte 

161  - 

240 

Text  (Table  A-1,  third  line) 

Note:  Latitude  of  center  (degrees) 
in  bytes  172-179;  Longitude  of 
center  (degrees)  in  bytes  189-196 

byte 

241  - 

320 

Blanks 

byte 

321  - 

400 

Text 

byte 

401  - 

480 

Blanks 

byte 

481  - 

560 

specification  of  first  data  channel 

(see  2.2) 


byte  3521  -  3600  specification  of  39th  data  channel 


2.2  Contents  of  specification  field 


Each  specification  field  mentioned  above  is  80  bytes  in 
length,  and  describes  one  data  channel.  The  format  is 
currently  as  follows: 


byte  2-3 
byte  7-8  : 

byte  13-15: 
byte  16-23: 

byte  24-31: 

byte  32-40: 
byte  41-60: 
byte  61-80: 


channel  Identification  (e.g.,  21) 
site  name  (e.g.,  FO) 
component  (e.g.,  SPZ) 
position  relative  to  center  of  array 
(NS),  in  meters 

position  relative  to  center  of  array 
(ES),  in  meters 

elevation  (meters  above  sea  level) 
Instrument  type  (e.g.,  GS-13) 
comments 


The  ARRAY  GEOMETRY  RECORD  is  distributed  on  all  tapes  in 
either  ASCII  or  EBCDIC.  The  record  gives  the  array  configura¬ 
tion  corresponding  to  the  actual  data  on  tape.  An  example 
is  given  in  Table  A-1. 


3.  DATA: 


Each  one  second  data  block  Is  3540  bytes  In  length  and 
consists  of: 

RECORD  HEADER  40  bytes 

1  status  frame  100  bytes 

34  data  frames  34*100  bytes 

RECORD  HEADER: 

byte  1  -  2  EBCDIC:  data  identification  EBCDIC  'NO', 

(CFDO)  Hexadecimal  X'CFDO' 

byte  3  -  4  Binary:  sampling  delay  (LP)  in  tenths  of  ms 

(currently  1008,  i.e.,  100.8  ms) 
byte  5-6  Binary:  sampling  delay  (IP)  in  tenths  of  ms 

(currently  508,  i.e.,  50.8  ms) 
byte  7-8  Binary:  sampling  delay  (LP)  in  tenths  of  ms 

(currently  60,  i.e.,  6.0  ms) 
byte  9-10  Binary:  correction  in  milliseconds  (fixed) 

for  UTC.  Current  value  is  1000,  i.e., 
1000  ms  should  be  suotracted  from 
all  data  times 
byte  11-40  Binary:  spare 

DATA: 

35  frames  (1  status  frame,  34  data  frames)  as  described  in 
HUB  documentation.  (Table  A-2) 

STATUS  FRAME:  word  50  of  this  frame  contains  the  total 
number  of  frames  in  this  record  (1-35). 

Note  that  many  conventions  used  for  NORESS  status  Information 
and  data  follow  those  established  for  the  RSTN  network. 
Reference  is  made  to  Breding  (1983)  for  additional  detail. 


Reference: 

Breding,  D.  (1983);  Data  users'  guide  for  the  Regional 

Seismic  Test  Network  (RSTN),  Sandla  report  SAND82-2935. 
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Table  A-1 

Example  of  an  Array  Geometry  Record. 
Each  line  corresponds  to  80  bytes. 


AG  - 

ARRAY 

GEOMETRY  RECORD 

09/27/85 

13:30  GMT 

NO  - 

NORESS 

DATA 

40. 

.OHZ  , 

10. OHZ  AND  . 

l.OHZ 

FO  - 

LAT(D) 

60.7353  LON(D)  11. 

5414 

CID 

SITE 

COMP 

NS(M) 

EW(M) 

ELEV(M) 

INSTRUMENT 

21 

FO 

SPZ 

0 

0 

247 

S-3 

61 

FO 

SPN 

0 

0 

247 

S-3 

A1 

FO 

SPE 

0 

0 

247 

S-3 

22 

AO 

SPZ 

3 

4 

302 

GS-13 

TEMP, 

62 

AO 

SPN 

3 

4 

302 

GS-13 

TEMP 

A2 

AO 

SPE 

3 

4 

302 

GS-13 

TEMP 

OD 

A1 

SPZ 

146 

49 

291 

GS-13 

OE 

A2 

SPZ 

-  103 

108 

311 

GS-13 

17 

B1 

SPZ 

321 

70 

299 

GS-13 

18 

B2 

SPZ 

30 

334 

315 

GS-13 

OA 

B3 

SPZ 

-  298 

143 

314 

GS-13 

OB 

B4 

SPZ 

-  217 

-  228 

299 

GS-13 

OC 

B5 

SPZ 

163 

-  272 

289 

GS-13 

13 

Cl 

SPZ 

687 

109 

299 

GS-13 

OF 

C2 

SPZ 

341 

603 

339 

GS-13 

TEMP 

14 

C3 

SPZ 

-  238 

647 

352 

GS-13 

23 

C4 

SPZ 

-  657 

208 

311 

GS-13 

63 

C4 

SPN 

-  657 

208 

311 

GS-13 

A3 

C4 

SPE 

-  657 

208 

311 

GS-13 

15 

C5 

SPZ 

-  569 

-  396 

299 

GS-13 

16 

C6 

SPZ 

-  48 

-  687 

303 

GS-13 

24 

C7 

SPZ 

548 

-  447 

275 

GS-13 

64 

C7 

SPN 

548 

-  447 

275 

GS-13 

A4 

C7 

SPE 

548 

-  447 

275 

GS-13 

04 

D1 

SPZ 

1480 

192 

305 

GS-13 

05 

D2 

SPZ 

1015 

1098 

372 

GS-13 

06 

D3 

SPZ 

76 

1493 

453 

GS-13 

07 

D4 

SPZ 

-  901 

1189 

379 

GS-13 

08 

D5 

SPZ 

-1451 

335 

348 

GS-13 

09 

D6 

SPZ 

-1326 

-  681 

352 

GS-13 

10 

D7 

SPZ 

-  566 

-1368 

337 

GS-13 

11 

D8 

SPZ 

414 

-1336 

301 

GS-13 

12 

D9 

SPZ 

1257 

-  802 

278 

GS-13 

Cl 

EO 

IPZ 

2 

3 

247 

KS-36000-04 

Cl 

EO 

IPN 

2 

3 

247 

KS-36000-04 

Cl 

EO 

IPE 

2 

3 

247 

KS-36000-04 

Cl 

EO 

LPZ 

2 

3 

247 

KS-36000-04 

Cl 

EO 

LPN 

2 

3 

247 

KS-36000-04 

Cl 

EO 

LPE 

2 

3 

247 

KS-36000-04 

COMMENTS 


POSITION 

POSITION 

POSITION 


CID 
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Table  A- 2 

HUB/Element  Data  Link.  Interface 


HUB  is  equipped  with  26  synchronous  2400  bit/sec  ports  for  element 
frame  acquisition.  Each  element  uses  the  HUB  transmitter  (element 
receiver)  clock  as  its  data  frame  transmitter  clock.  At  initialization 
each  element  will  enable  output  starting  when  the  last  byte  of  a  com¬ 
mand  frame  is  received  from  HUB.  This  method  of  data  frame  synchroni¬ 
zation,  which  applies  to  all  elements,  defines  a  single  sampling 
interval  throughout  the  array,  so  time  correlation  of  individual  ele¬ 
ment  data  is  simple.  RSAS  elements  generate  50  word  Short  Period  (SP) 
or  Broad  Band  (BB)  frames.  Three  axis  elements  generate  3  SP  frames/ 
sec,  one  for  each  axis.  Formats  and  data  definitions  for  the  2  frame 
types  are  given  in  the  article  below. 

SP  Frame  Format 


15|  14|  13|  12|  111  10|  09|  08|  07|  06|  05|  04|  03|  02|  01|  00| 


FRAME  SYNC  WORD 

1  (F325  HEX) 

FRAME  SYNC  WORD 

2  (OCDA  HEX) 

0  1 

2  1 

TIME  WORD  1 

1  15| 
f2  1 

16| 

2  1 

TIME  WORD  1  1  23|  07|  CHANNEL 

|2  |2  1 

STATUS  (CS) 

1  o| 
|2  1 

40  WORDS  OF  SHORT  PERIOD  DATA 

CHANNEL  SOH 

13| 

2  1 

COMMAND 

1  o| 
|2  1 

COMMAND 

1  AID 

161  03|  1 

2  |2  I  |2 

1  TIME  WORD  2 

oj  27 1  26 1  25 1 

|2  |2  |2  (2 

24|  19|  18| 

|2  |2  |2 

1 

17|  16| 

|2  1 

TIME  WORD  2 

1  oo| 
|2  1 

CID 

1  FA 

HEX 

6  BITS 

UNDEFINED  (*)  | 

10  BIT  AUTHENTICATION  NUMBER 

*  The  upper  6  bits  of  any  authentication  word  must  be  set  to  0 
before  authentication  comparison  is  performed. 


SP  frame  channel  status  format 


Each  Z ,  E  or  N  SP  frame  contains  one  Channel  Status  (CS)  byte  con¬ 
taining  digital  status  for  the  axis  as  defined  in  the  table  below. 

Bit  Function 

0  unused 

1  CAL  flag 

2  Clip  flag 

3  mass  motor  busy 

4  period  motor  busy 

5  mass  position  motor  busy 

6  motor  direction 

7  unused 


SP,  BB  frame  channel,  frame  ID  (CID)  format 

Each  SP  or  BB  frame  CID  byte  contains  a  6-bit  channel  and  a  2-bit 
frame  type  identification  field  as  shown  below.  The  CID  byte  is  a 
constant  programmed  by  a  strap  port  in  each  element.  For  SP  and  BB 
elements  only  1  CID  byte  is  required.  A  TSP  element  generates  3  CID 
bytes  and  the  strap  value  will  be  assigned  to  the  Z  axis  SP  frame.  The 
element  processor  will  add  40  hex  or  80  hex  to  the  Z  axis  CHANNEL  ID 
for  the  N  and  E  axis  SP  frame  CID  byte,  respectively. 


SP,  BB  frame  time  word  formats 

A  28-bit  frame  count  incremented  once  per  second  and  used  for  authen¬ 
tication  variable  selection  by  the  element  processor  is  the  source  of 
TIME  WORD  1  and  2  for  all  SP  and  BB  frames.  The  Integrity  of  the  24 
most  significant  bits  of  frame  count  is  supported  by  three  8-bit  ports 
which  have  battery  backup  for  power.  TIME  WORD  1  contains  bits  2**0 
through  2**23  of  frame  count  in  reverse  order.  TIME  WORD  2  contains 
frame  count  bits  2**27  through  2**24  and  bits  2**19  through  2**0. 


Seismic  data  format 

The  least  significant  14  bits  of  any  seismic  data  word  are  the  digi¬ 
tized  representation  of  the  sensor  analog  data  signal  and  the  1  most 
significant  bits  of  the  word  are  encoded  by  the  element  processor  to 
indicate  the  GAIN  SELECT  value  as  follows: 

BIT  15  14 

0  0  GAIN  SELECT  =  128 
0  1  GAIN  SELECT  -  32 
1  0  GAIN  SELECT  =  8 
1  1  GAIN  SELECT  -  1 


FRAME  SYNC  WORD  1  (F325  HEX) 


FRAME  SYNC  WORD  2  (OCDA  HEX) 


0|  TIME  WORD  1  f  15| 

2  |2  I 


16 1  TIME  WORD  1  |  23|  07|  BB  CHANNEL  STATUS  |  0| 

2  I  |2  |2  I  |2  I 


10  WORDS  OF  IPZ  DATA 
APPEARING  IN  ORDER  COLLECTED 


10  WORDS  OF  IPN  DATA 
APPEARING  IN  ORDER  COLLECTED 


10  WORDS  OF  IPE  DATA 
APPEARING  IN  ORDER  COLLECTED 


LPZ  DATA  WORD 


LPN  DATA  WORD 


LPE  DATA  WORD 


GTSOH  STATUS 


DH  SOH  STATUS  WORD 


DH  SOH  STATUS  WORD 


DH  SOH  STATUS  WORD 


4  UNUSED  WORDS,  ALL  =  0 


15 1  COMMAND  |  0| 


2  1 

12  1 

COMMAND 

1  AID 

1  TIME  WORD  2 

19| 

161  03| 

1  oj  27)  26)  23)  24)  19)  18)  17)  16) 

2 

2  |2  1 

(2  |2  12  |2  |2  (2  {2  (2  |2  j 

15| 

TIME  WORD  2 

1  00) 

2  1 

|2  i 

CID 


FA  HEX 


6  BITS  UNDEFINED 


I  10  BIT  AUTHENTICATION  NUMBER 


HUB  status  frame  format 


i5| 

14|  13|  12|  11) 

lUj  09|  08|  07|  06|  05|  04|  03|  02|  Oij  00| 

FRAME 

SYNC 

WORD  1 

=  F325 

HEX 

FRAME 

SYNC 

WORD  2 

=  OCDA  HEX 

01  0| 

MASTER 

FRAME  NUMBER 

GMT  WORDS  1 

.  2.  3, 

4.  5 

HUB  ANALOG 

STATUS 

WORD  1 

THRU  10 

HUB  DIGITAL 

STATUS 

WORD 

PORT 

1  STATUS  BYTE 

1 

PORT 

0  STATUS  BYTE 

PORT 

3  STATUS  BYTE 

1 

PORT 

2  STATUS  BYTE 

PORT 

5  STATUS  BYTE 

1 

PORT 

4  STATUS  BYE 

PORT 

7  STATUS  BYTE 

1 

PORT 

6  STATUS  BYTE 

PORT 

9  STATUS  BYTE 

1 

PORT 

8  STATUS  BYTE 

PORT 

11  STATUS  BYTE 

1 

PORT 

10  STATUS  BYTE 

PORT 

13  STATUS  BYTE 

1 

PORT 

12  STATUS  BYTE 

PORT 

15  STATUS  BYTE 

1 

PORT 

14  STATUS  BYTE 

PORT 

17  STATUS  BYTE 

1 

PORT 

16  STATUS  BYTE 

PORT 

19  STATUS  BYTE 

1 

PORT 

18  STATUS  BYTE 

PORT 

21  STATUS  BYTE 

1 

PORT 

20  STATUS  BYTE 

PORT 

23  STATUS  BYTE 

1 

PORT 

22  STATUS  BYTE 

PORT 

25  STATUS  BYTE 

1 

PORT 

24  STATUS  BYTE 

GENERAL  STATUS  WORD 

UPLINK 

TIME 

POSITION  WORD 

LAST  HUB  COMMAND  WORD 

LAST 

COMMAND  BYTE 

1 

LAST 

COMMAND  ELEMENT 

13  WORDS  OF  MULTIPLEXED  DIGITAL  STATUS 

1  WORD,  CONTENT  UNDEFINED 

THIS  WORD  IS  RESERVED  FOR  SCARS  USE  AND  CANNOT  BE  USED  BY  HUB 


GMT  TIME  CODE  INTERFACE 


This  interface  consists  of  three  16-bit  ports  which  provide  time  of 
year  (GMT)  to  the  millisecond  which  is  inserted  into  HUB  status  frame 
GMT  words  1  through  5.  Format  of  each  word  is  shown  in  the  following. 


GMT  word  formats 


Bit 


Word  1 


Word  2 


Word  3 


0 

Units  of  miliisec. 

1 

Units  of  sec. 

1 

Units 

of  min. 

1 

1 

Units  of  miliisec. 

2 

Units  of  sec. 

2 

Units 

of  min. 

2 

2 

Units  of  miliisec. 

4 

Units  of  sec. 

4 

Units 

of  min. 

4 

3 

Units  of  miliisec. 

8 

Units  of  sec. 

8 

Units 

of  min. 

8 

4 

Tens  of  miliisec. 

1 

Tens  of  sec. 

1 

Tens 

of  min. 

1 

5 

Tens  of  miliisec. 

2 

Tens  of  sec. 

2 

Tens 

of  min. 

2 

6 

Tens  of  miliisec. 

4 

Tens  of  sec. 

4 

Tens 

of  min. 

4 

7 

Tens  of  miliisec. 

8 

0 

0 

8 

Hund.  of  miliisec. 

1 

0 

0 

9 

Hund.  of  miliisec. 

2 

0 

0 

10 

Hund.  of  miliisec. 

4 

0 

0 

11 

Hund.  of  miliisec. 

8 

0 

0 

12 

0 

0 

0 

13 

0 

0 

0 

14 

0 

0 

0 

15 

0 

0 

0 

iit 

Word  4 

Word  5 

0 

Units  of  hours 

1 

Units  of  days 

1 

1 

Units  of  hours 

2 

Units  of  days 

2 

2 

Units  of  hours 

4 

Units  of  days 

4 

3 

Units  of  hours 

8 

Units  of  days 

8 

4 

Tens  of  hours 

1 

Tens  of  days 

1 

5 

Tens  of  hours 

2 

Tens  of  days 

2 

6 

0 

Tens  of  days 

4 

7 

0 

Tens  of  days 

8 

8 

0 

Hund.  of  days 

1 

9 

0 

Hund.  of  days 

2 

10 

0 

0 

11 

0 

0 

12 

0 

±  1  millisecond  accuracy 

13 

0 

±  5  millisecond  accuracy 

14 

0 

±  5U  millisecond 

accuracy 

15 

0 

±  500  millisecond 

accuracy 

